1 


um 


Interim  Report 
CGC-RN105-77 


SUMMARY  AND  EXPERIMENTAL  PLAN 
FOR  REMOTE  SENSING  OF 
SUBSURFACE  TEMPERATURE  AND  SALINITY 


D.  A.  Leonard  and  B.  Caputo 

Computer  Genetics  Corporation 
Wakefield,  Massachusetts 


July  1977 


Prepared  F or: 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


niSTRIBUTTON  STATEMEWT  K_ 
Approved  tot  public  ratooa* 
Distribution  Unlinitid 


D D C 


APR  3 1979 


9 03  30  0.3 


T 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh#«  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


ACCESSION  N< 


^ND  EXPERIMENTAL  ,PLAN  FOR 

iMOTE  SENSING  OF  SUBSURFACE 

?EM  PE  RAT  U RE  AND  SALINITY-  

t _ a * — 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOG  NUMBER 


. tipeot  HEJSgi.A'gffliaB  eflvgHPP 


AUTHOR!  «J 


D.  A. /Leonard  — M B.^apato  J 


PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Computer  Genetics  Corporation  / 
18  Lakeside  Office  Park 
Wakefield,  Massachusetts  01880 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  Virginia  ZZ217 


1*  MONITORING  AGENCY  NAME  ft  AODRESSfl*  dlllormnt  /rom  Controlling  Ollleo) 


Interim  ^epGft^e 

1 Jul4p&31  Oct«*76 
PeUPWJi'Wgllli.'HBijSTwuwi 


6.  CONTRACT  OR  CR*NT  NUMBERf 

ijXfjsja  14- 76 -C- IjSfs'l 1 (-4 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  & WORK  UNIT  NIIMRERS 


NUMBER  OF  PAGES 

105 


IS.  SECURITY  CLASS,  (a!  IM,  r«.portJ 

UNCLASSIFIED 


IS«.  DECLASSIFICATION'  DOWNGRADING 
SCHEDULE  N/A 


16  DISTRIBUTION  STATEMENT  fo!  ll.lt  R*port> 


PIBTIUBUTION  STATEMPtT  K 
Approved  foi  public  isleoiBj 
Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  the  •batrmct  entered  In  Block  70,  If  dllloronl  I rom  Report) 


Ift.  SUPPLEMENTARY  notes 


19  KEY  WORDS  (Contlnuo  on  / eeere e tide  It  noceeeary  and  Identity  try  block  numbtf) 

Water  Temperature 
Salinity 
Subsurface 
Remote  Sensing 
Raman  Lidar 


l^ft^VeSTR  AcT (Contlnuo  on  reveree  rid*  II  necoeeary  and  Idontlly  by  block  nuaibar) 

The  temperature  and  salinity  of  liquid  water  can  be  inferred  from  its 
Raman  spectrum.  This  report  assesses  and  evaluates  the  prominent  factors) 
which  influence  the  ability  of  a remote  sensing  Raman  lidar  to  measure  sub- 
surface profiles  of  water  temperature  and  salinity.  A plan  is  developed  for 
the  experimental  investigation  of  the  feasibility  of  the  technique  in  ocean  and 
coastal  waters.  Interferences,  including  surface  wave  effects,  differential 
attenuation,  differential  depolarization  and  fluorescence  are  considered.  * — 
over  I 


DD  , 


FORM 
JAN  7) 


1473 


EDITION  OF  I NOV  *S  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (*h«o  Dm,  E nt,r*d) 


3/J.  r 73 


TABLE  OF  CONTENTS 


Section  Page 

1.0  INTRODUCTION  AND  SUMMARY  1 

2.0  PHYSICAL  BASIS  OF  EXPERIMENT  AND  PRIOR  WORK  2 

3.0  NEW  WORK  ACCOMPLISHED  UNDER  SUBJECT  CONTRACT  5 

3.  1 Major  Interferences  and  Noise  Sources  5 

3.  1.  1 Surface  Wave  Effects  5 

3.  1.1.1  Differential  Transmission  of  Orthogonal  6 

Linear  Polarization 

3.  1.1. 2 Depth  Resolution  12 

3.  1.1. 3 Lateral  Resolution  14 

3.1.2  Volume  Transmission  Phenomena  17 

3.  1.2.1  Differ ential  Depolarization- Experimental  17 

3. 1.2. 2 Differ  ential  Depolarization-Theoretical  21 

3.  1.2.3  Differential  Attenuation  and  the  Need  for  24 

Multiwavelength  Polarization  Detection 

3.1.3  Non-Raman  Photon  Sources  28 

3.  1.3.  1 Ambient  Light  Background  28 

3.  1.3.  2 Laser  Induced  Fluorescence  29 

3.2  Detailed  Experimental  Apparatus  Design  30 

3.2.1  Laser  Raman  Spectrometer  for  Field  30 

Utilization 

3.2.2  Research  Vessel  Facility  32 

3.2.3  Preliminary  Experimental  Raman  Data  37 

Obtained  From  Natural  Waters 

3.3  Experimental  F ield  T est  Plan  41 

3.3.  1 Experimental  Objectives  and  Methods  43 

3.  3.  1.1  Transmis sion  (i.  e.  K)  43 

3.  3.  1.2  Fluorescence  to  Raman  Ratio  (F)  43 

3.3.  1.3  Differential  Polarization  Analysis  43 

3.3. 1.4  Proportionality  Between  and  K 43 

3.3. 1.5  Temperature  vs.  Salinity  Separability  44 

3.3.2  Sea  Truth  Requirements  44 

3.3.2.  1 Sea  Truth  Measurement  Parameters  45 

3.3.3  Test  Site  Requirements  46 

3.3.3.  1 Water  Quality  46 

3.3.  3.2  Salinity  47 

3.4  Q & A on  Fundamental  and  Planned  Experimental  62 

Techniques 


i 


table  of  contents 

(Continued) 


Section  Page 

4.0  CONCLUSIONS  AND  RECOMMENDATIONS  67 

REFERENCES  70 

APPENDIX  A - Physical  and  Mathematical  Justification  A1  - A8 

APPENDIX  B - Computer  Program  to  Determine  the  B1  - B5 

Effect  of  Ocean  Surface  Wave  Shape  on 
Differential  Transmission  of  Orthogonal 
Polarizations 

APPENDIX  C - Computer  Program  to  Determine  the  Cl  - C8 

Effect  of  Ocean  Surface  Waveshape  on 
Refractive  Depth  Penetration  Correction 

APPENDIX  D - Depolarization  Calculations  D1  - D3 

APPENDIX  E - Experimental  Remote  Sensing  of  El  - E4 

Subsurface  Temperature  In  Natural 
Ocean  Water  (CGC  Research  Note  102) 


LIST  OF  ILLUSTRATIONS 


Figure 


3-14 


3-18 

3-19 

3-20 


3-22 


Air-Water  Interface  Geometry  in  the  Plane  of  Incidence 

Differential  Polarization  vs.  Incident  Angle 

Differential  Transmission  of  Orthogonal  Linear 
Polarizations  Averaged  Over  Various  Wave  Spectra 

Geometrical  Optics  for  Wave  Refraction  Depth  Correction 
and  Lateral  Error  Correction 

Fractional  Error  in  Depth  Measurement  Due  to  Refraction 
as  Function  of  H/L  for  Stokes,  Trochoidal  and  Sinusoidal 
W aves 

Fractional  Error  in  Depth  Measurement  Due  to  Refraction 
as  Function  of  H/L  for  Stokes  Waves  Showing  Average, 
90%  and  Maximum  Values 

Depolarization  During  Transmission  Through  Natural 
Waters  Obtained  by  S.  Q.  Duntley 

Differential  Attenuation  vs.  Wavelength  Separation 
With  Depth  as  a Parameter  (K  Minimum) 

Differential  Attenuation  vs.  Wavelength  Separation 

A Frc 


With  Depth  as  a Parameter  (100  A From  K Minimum) 

Functional  Block  Diagram  of  Dual  Polarization  Raman 
Spectrometer 

Bistatic  Laser  Raman  Equipment  Aboard  the  Research 
Vessel  Makai 

Block  Diagram  of  Laser  Raman  System  Currently 
Operational  on  Research  Vessel  Makai 

Bistatic  Laser  Raman  Equipment  in  the  Laboratory 

Bistatic  Laser  Raman  Equipment  Aboard  the  Research 
Vessel  Makai 

Raman  Spectral  Data  for  Liquid  Water 

Schedule  and  Plan  for  Raman  Temperature /Salinity  Research 

Map  of  Northern  Massachusetts  and  New  Hampshire 
Coastal  Region 

R/V  Berth  Near  Blyman  Canal 

R/V  Berth  and  Annisquam  River 

Gloucester  Harbor,  Raw  Sewage  Outlet 

Deep  (56 ' ) "Clear"  Water  Hole  Off  Coffins  Beach 

Deep  Open  Waters  Off  Castle  Neck 

Merrimack  River  Test  Site 


3-24  Sheltered  Shakedown  Area 


A-l  Raman  Spectrum  and  Structure  of  Liquid  Water  A2 


1 Block  Diagram  of  Laser  Raman  System  Currently  E3 

Operational  on  Research  Vessel  Makai 

2 Raman  Spectral  Data  for  Liquid  Water  E3 

3 Typical  System  Operating  Map  (100  Meter  Mixed  Layer)  E4 


LIST  OF  TABLES 

TabIe  Page 

2- 1  Comparison  of  Raman  and  Volume  Scattering 

Coefficients  for  Water 

3- 1  Sample  Calculation  Results  for  Depolarization 

During  Transmission 

3-2  Water  Discharge  Data  of  Merrimack  River  Below 

Concord  River  at  Lowell,  Massachusetts  in 
CFS  (cubic  feet  per  second) 

3-3  Climatological  Table 

3-4  Mean  Surface  Water  Temperatures  (T)  and  Densities  (D) 

3-5  Meteorological  Table  for  Coastal  Area  Off  Boston 

A-l  Baseline  System  Parameters 

A-2  Underwater  Windows  for  Various  Ocean  Waters 


rrr — ” — - : — 

1.0  INTRODUCTION  AND  SUMMARY 

• This  report  summarizes  the  work  accomplished  on  Contract  N00014-76-C- 

1007  Task  I during  the  period  1 July  1976  to  31  October  1976.  The  effort 
concentrated  on  planning  and  studying  the  experimental  investigation  of  the 
feasibility  of  applying  Raman  spectroscopy  to  obtain  subsurface  temperature 
and  salinity  information  in  ocean  and  coastal  waters. 

A critical  review  and  study  was  made  of  the  proposed  experimental  plan 
described  in  Computer  Genetics  Corporation  proposal  CGC-P121-76  dated 
April  1976.  The  calculations  and  assumptions  upon  which  the  experiment 
is  based  were  expanded,  reviewed  and  with  the  exception  of  a newly  measured 
and  larger  value  for  the  basic  cross-section,  were  found  to  be  consistent 
with  current  information. 

Much  of  the  new  work  performed  on  Contract  N00014-76-C-1007  centered 
on  the  identification  and  quantification  of  the  various  interferences  and  noise 
sources  that  could  limit  system  performance  and  which  have  been  quantitatively 
evaluated  with  reassuring  results. 

Areas  requiring  new  and  supplemental  experimental  attention  to  develop 
the  technique  have  been  identified  and  a broad  as  well  as  detailed  experimental 
plan  has  been  written.  An  existing  government  owned  apparatus  configuration 
was  used  to  obtain  the  first  Raman  data  spectrum  taken  from  a boat  in  the 
natural  coastal  environment.  The  temperature  obtained  from  the  spectrum 


agreed  with  prior  laboratory  calibration  and  the  measured  fluorescence  to 
Raman  ratio  was  very  much  less  than  unity.  A preferred  configuration  which 
can  be  readily  assembled  from  additional  existing  government  inventory  is 
presented. 

The  refractive  effects  of  surface  waves  related  to  differential 
polarization  transmission  and  depth  correction  were  studied  and  computer 
programs  were  written  and  exercised  which  show  less  than  a 1.0  C offset 
in  the  temperature  measurement  and  less  than  a few  ppt  in  salinity 
measurements  even  if  uncorrected.  The  major  limiting  factors  in  depth 
resolution  have  been  identified  as  laser  pulse  duration  and  the  random 
refractive  effects  at  the  air  water  interface.  The  former  effect  is  reduced 
with  advanced  laser  design,  while  through  computer  calculation  the  latter 
effect  appears  to  stabilize  as  a 0.  1 to  0.3  percent  effect. 
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Depolarization  during  transmission  through  natural  waters  has  been 
studied  and  the  concensus  based  upon  prior  work  of  several  investigators 
is  that  over  80%  of  the  original  polarization  remains  after  traversing  20 
attenuation  lengths.  Also  it  has  been  calculated  that  if  the  proportionality 
between  the  depolarization  coefficient  p g and  the  diffuse  attenuation  coefficient 
K can  be  known  to  within  10%,  that  after  corrections  are  applied,  accuracies 
of  0.  16°C  in  temperature  and  0.  7 ppt  in  salinity  per  attenuation  length  can 
be  expected.  This  means,  for  example,  that  at  a depth  of  three  attenuation 
lengths  a temperature  uncertainty  of  0.  48°C  and  a salinity  uncertainty  of 
2.  1 parts  per  thousand  would  be  expected.  The  range  of  temperature  and 
salinity  considered  was  - 2°C  to  35°C  and  0 ppt  to  40  ppt  respectively. 

Through  computer  simulation,  tables  have  been  generated  to  relate  the 
theoretical  accuracy  of  the  polarization  ratio  measurement  to  spectrometer 
bandwidth  and  multiwavelength  detection  as  a function  of  differential 
attenuation. 

A well  researched  and  documented  initial  test  site  which  offers  a 
representative  cross-section  of  water  quality  and  natural  coastal  effects 
such  as  sediment  loading,  salt  water  intrusion,  thermal  gradients,  etc.  is 
discussed  in  detail. 

This  Technical  Summary  Report  highlights  the  results  of  the  effort. 
Supporting  technical  references,  computer  calculation  printouts,  design 
information,  component  specification  and  cost  backup  is  on  file  in  various 
formats  at  the  CGC  office  in  Wakefield,  Massachusetts. 

2.  0 PHYSICAL  BASIS  OF  EXPERIMENT  AND  PRIOR  WORK 

The  phy sical  basis  for  the  Raman  measurement  of  subsurface  temperature 
and  the  prior  work  upon  which  the  current  program  is  based  were  described 
in  Computer  Genetics  Proposal  CGC-P121-76  dated  April  1976.  Appendix  A 
of  this  report  repeats  the  basic  information  relating  to  the  physical  basis 
and  mathematical  justification  as  it  appeared  in  that  proposal. 

During  the  course  of  the  literature  search  performed  as  part  of  the 

new  work  conducted  on  Contract  N00014-76-C -1007  another  more  recently 

measured  value  for  the  basic  liquid  water  Raman  cross-section  was 

discovered.  A paper  entitled,  "Temperature  Dependence  and  Cross- 

Sections  of  Some  Stokes  and  Anti-Stokes  Raman  Lines  in  Ice"  by  Slusher 

(and  Derr^  also  reports  an  experimentally  measured  value  for  the  liquid 

-33  2 

water  Raman  cross-section  of  4.  5 x 10  m /molecule  sr.  This  cross- 
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section  value  is  approximately  a factor  of  five  higher  than  the  value  that 
appeared  in  an  older  reference'  ' and  which  has  previously  served  as  the 
basis  for  our  system  performance  calculations. 

In  a private  conversation  with  Derr  concerning  the  liquid  water  cross- 

section  measurements  he  expressed  the  opinion  that  he  was  not  aware  of 

(2) 

the  prior  work.  He  did  feel  that  his  own  work  was  accurate  to  probably 
much  better  than  a factor  of  five.  Derr  also  pointed  out  that  it  is  common 
for  liquid  phase  Raman  cross-sections  to  be  much  larger  than  gas  phase  even 
on  a per  molecule  basis  and  the  fact  that  this  is  true  for  water  was  not 
surprising. 

Table  2-1  shows  a representative,  but  by  no  means  complete,  compilation 
of  volume  backscattering  coefficients  that  have  been  measured  by  various 
investigators  for  distilled  water  and  ocean  and  lake  water  of  various  types. 

Also  included  in  Table  2-1  are  Raman  backscattering  coefficients  for  distilled 
water  and  ice.  The  wavelength  of  light  used  for  the  measurements  cited  in 
Table  2-1  varies  somewhat  from  investigator  to  investigator,  but  all  use 
light  in  the  blue-green  near  500  nm  so  that  order  of  magnitude  comparisons 
may  be  considered  valid. 

From  the  data  compiled  in  Table  2-1  it  can  be  noted  that  the  new  Slusher 
and  Derr  value  for  the  Raman  coefficient  is  essentially  equal  to  the  on-frequency 
backscatter  coefficient  of  Tyler  for  distilled  water.  Far  from  being  a "weak" 
effect  the  Raman  is  comparable  to  the  main  on-frequency  return  from  the 
water  itself. 

This  has  implications  when  making  comparisons  between  the  Raman 
method  and  other  laser  scattering  methods  for  remote  subsurface 
measurements  such  as  the  recently  suggested'  Brillouin  measurements 
to  determine  sound  speed. 

It  should  also  be  noted  that  the  cross-section  fo-  ice  is  comparable  to 
that  of  liquid  water.  Slusher  and  Derr^1*  have  obtained  temperature 
measurements  using  Stokes  to  anti-Stokes  line  ratios  in  ice  from  0°C  to 
-50°C.  A Raman  instrument  capable  of  measuring  temperature  and  salinity 
should  also  be  capable  not  only  of  detecting  ice  but  also  of  measuring  ice 
tempe  ratures. 


Comparison  of  Raman  and  Volume  Scattering  Coefficients  for  Water 


Raman  Scattering  Coefficient,  (X=  488  nm) 

(Meters  * Steradian  *) 

Distilled  Water  . 00003  Chang  and  ~ 

Young(  1 974) 

.00015  Slusher  and 

De  rr(  1 975)  1 

Ice  . 000093  Slusher  and 

De  r r ( 1975)  1 


Volume  Scattering  Coefficient,  a (180°) 
(Meters  * Steradian  ) 


Distilled  Water 

. 0001 7-. 00020 

Tyler(l  96 1 ) 3 

Ocean 

Various  Locations,  Eastern  North 
Atlantic 

. 0005 

J e rlov(  1 96 1 )4 

Off  Bermuda  (31°  57'N,  65°  11'  W) 

. 002 

Morrison(1970)^ 

Coastal 

Between  Woods  Hole  and  the  180  m 
Contour 

.003-. 004 

Spilhaus(l  968)  ^ 

Near  Santa  Catalina  Island 
(33°  32'N,  118°  17'W) 

. 00020 

Tyler(l  9 6 1 )3 

Long  Island  Sound  (41°  16'N,  70°  08'W) 

. 02 

Morrison(  1 970)3 

Lake 

Winnipesaukee,  New  Hampshire 

Pend  Oreille,  Idaho 

. 003 

. 00 1 7- . 0041 

Duntley(  1 963)^ 
Tyler(196l)3 
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3.0  NEW  WORK  ACCOMPLISHED  UNDER  SUBJECT  CONTRACT 


This  section  of  this  Technical  Summary  Report  summarizes  the  work 
accomplished  on  Contract  No.  N00014-76-C-1007  during  the  four  month 
period  1 July  1976  to  31  October  1976.  The  effort  concentrated  on  planning 
and  studying  the  experimental  investigation  of  the  feasibility  of  applying 
Raman  spectroscopy  to  obtain  subsurface  temperature  and  salinity  information 
in  natural  ocean  and  coastal  waters. 

The  effort  included  literature  searches,  interviews  with  key  personnel 
in  the  ocean  sciences  area,  theoretical  calculations  including  the  development 
and  utilization  of  several  new  computer  routines,,  a survey  of  the  availability 
and  the  coordination  of  the  utilization  of  Government  property  in  a laser  Raman 
coastal  experiment,  the  planning  of  specific  experimental  apparatus  designs 
and  specific  test  plans,  and  the  analysis  and  interpretation  of  preliminary 
field  data. 

This  1 echnical  Summary  Report  highlights  the  results  of  the  effort. 
Supporting  technical  references,  computer  calculation  printouts,  design 
information,  component  specification  and  cost  backup  are  on  file  in  various 
formats  at  the  CGC  office  in  Wakefield,  Massachusetts. 

3.  1 Major  Interferences  and  Noise  Sources,  An  important  aspect  of 
the  effort  was  to  identify  and  quantify  the  major  interferences  and  noise 
sources  that  laser  Raman  spectroscopic  measurements  of  subsurface  ocean 
temperature  and  salinity  would  encounter.  The  purpose  was  to  ensure  that 
subsequent  experiments  would  be  properly  designed  and  planned  so  that  the 
significant  questions  concerning  interferences  and  noise  sources  could  be 
adequately  answered. 

The  following  subsections  summarize  the  interferences  and  noise 
sources  as  follows:  Section  3.  1.  1 describes  surface  wave  effects  both  as  they 
influence  the  depth  resolution,  and  also  as  a depolarization  noise  source; 

Section  3.  1.  2 describes  the  volume  transmission  phenomena  of  differential 
attenuation  and  differential  depolarization;  and  Section  3.  1.  3 describes  the 
influence  of  non-Raman  photon  sources  such  as  fluorescence  and  ambient 
backgrounds. 

3.  1.  1 Surface  Wave  Effects.  The  effect  of  surface  waves  were  found 
to  relate  in  three  significant  ways  to  the  subsurface  measurement  of  tempei  ture 
and  salinity.  These  are  (a)  the  differential  transmission  of  orthogonal  linear 
polarizations  through  the  air-water  interface,  (b)  the  loss  of  depth  resolution 
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due  to  variable  refraction  angles  produced  at  a wavy  water  surface  and 
(c)  lateral  error  due  to  the  same  variable  refraction  angles. 

3 . 1 . 1.1  Differential  Transmission  of  Orthogonal  Linear  Polarization. 
The  differential  transmission  of  orthogonal  linear  polarizations  through  an 
oblique  air-water  interface  in  terms  of  typical  wave  spectra  was  investigated 
by  means  of  "exact"  one -dimensional  calculations  on  an  IBM  360  computer. 
New  computer  routines  were  written  which  used  the  basic  Fresnel  refraction 
equations.  ^ Sinusoidal  as  well  as  finite  amplitude^1  trochoidal  and 
"Stokes"  waves  were  employed  in  the  program.  The  effect  of  variable 
crest  height  to  wavelength  ratios  was  investigated  in  order  tp  study  the 
effect  of  sea  state  on  the  measurement. 

The  results  showed  that  even  at  the  worst  case  sea  state 
condition  of  a Stokes  limit^10*  wave  amplitude  peak-to-peak  to  wavelength 
(H/L)  ratio  of  0.  143  the  effect  is  to  introduce  a bias  of  less  than  1.0%  in  the 
Raman  depolarization  ratio  measurement.  This  would  mean  less  than  a 
1.0°C  off-set  in  the  temperature  measurement  or  a few  parts  per  thousand 
error  inthe  salinity  measurement  if  uncorrected.  It  was,  therefore, 
concluded  that  in  the  initial  experimental  configuration  no  attempt  should  be 
made  to  improve  the  accuracy  by  correlation  on  a pulse-to-pulse  basis  of 
the  wave  angle,  an  average  over  a given  wave  spectra  being  adequate  for 
initial  experimental  investigation. 


Figure  3-1.  Air-Water  Interface  Geometry  inthe  Plane 
of  Incidence. 
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The  details  of  the  calculation  are  illustrated  in  Figure  3-1  which  shows 
the  air-water  plane  of  incidence  with  the  incident  light  of  intensity  A shown 
at  an  angle  0^  with  respect  to  the  normal  to  the  air- water  interface.  The 
reflected  ray  R and  the  refracted  ray  B at  angle  0 are  also  shown. 

1 (9) 

Using  the  Fresnel  refraction  equations  in  the  form: 


tan  (©  - © ) 

3.  - i 5 

tan2  (©.  + ©J 
l t 


R±  = 


sin2  ( 0.  - ©t) 


sin  ( ©i  + ©fc) 


the  differential  transmission  of  polarization  through  the  surface,  i.  e.  B,  - 
/B^  = AB/Bj^  can  be  easily  derived  from  the  conservation  equations 
B = 1 - R and  B.  = 1 - R as  follows: 

//  O 1 1 


tan  ( © - © ) 

R . = 1 - R,  =1 = — 

/ X tan2  ( ©.  + © ) 

i t 


B = A - R = 1 - R = 1 

111  1 


sin2  ( ©.  - © ) 
l t 

sin  ( ©.  + ©.) 

l t 


AB  = ^ - Bx 


sin2  (©.  - © ) 
l t 

sin2  (©.  + © ) 
l t 


tan  (©.  - ©,) 
l t 

2 

tan  (©.  + ©J 
l t 


sin2  (0i  - 0t) 
sin2  (©j.  + 0t) 


tan  (B.  - 0t) 
tan2  (0.  + 0fc) 

sin2  (Qi  - 0t) 
sin2  (©i  +©t) 
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BJ_  sin2  (0.  +0t)  - sin2  (0i  -0t) 


The  effect  of  oblique  transmission  through  an  air-water  interface  is 
thus  to  produce  differential  transmission  of  orthogonal  linear  polarizations. 
If  the  incident  light  is  a pure  circular  polarization  mode,  the  effect  will  be  to 
produce  elliptical  light. 

The  quantitative  nature  of  the  effect  is  shown  in  Figure  3-2,  which  is  a 
plot  of  the  differential  polarization  produced  by  oblique  transmission  through 
an  air-water  interface,  vs.  the  incident  angle,  As  shown,  the  effect  is 

"negligible"  for  small  angles  and  becomes  nearly  unity  at  angles  approaching 


ial  Polarization  vs.  Incident  Angle. 


In  a practical  situation,  the  average  differential  polarization 
obtained  by  integration  over-typical  wave  spectra  must  be  considered. 
Calculations  were  made  using  sinusoidal  and  finite  amplitude  trochoidal 
and  Stokes  waves  of  varying  crest  height  to  wavelength  ratio.  The 
mathematical  formulation  of  the  trochoidal  and  Stokes  waves  was  taken 
from  Kinsman^). 

The  listing  of  the  above  described  computer  program  which  was 
developed  as  part  of  the  subject  contractual  effort  is  exhibited  as  Appendix  B 
of  this  report  and  is  given  the  name  SWPUN  (Surface  Wave  Polarization 
Unbalance).  Also  exhibited  in  Appendix  B are  typical  computer  printouts 
showing  the  results  of  calculations. 

The  calculational  procedure  for  the  sine  wave  case  is  shown  below 
by  way  of  example.  The  more  complex  trochoidal  and  Stokes  wave  calculations 
can  be  derived  from  the  listing  in  Appendix  B. 

y = A sin  2 jt  — ^ — 

lml 

At  any  point  on  the  wave  tan  = dy/dx  or 

©i  = arc  tan  ( 2jt  — — cos  2ir  ). 


Using  in  addition,  Snell's  Law  relation  sin  0.  = n sin  ©.  the  value  of 

X l 

AB/B^may  be  computed  for  any  value  of  A/L  over  the  entire  sine  wave. 

An  effective  AB/B^  averaged  over  the  sine  wave  is  thus  calculated  where 

100 

J = 1 

for  O < < 1 with  A (j)  = 0.01 

Figure  3-3  shows  the  results  of  the  calculations  displayed  as  a plot  of 
the  effective  differential  polarization  transmission  ( AB  /B±  )eff  as  a function  of 
H/L,  the  peak  to  trough  wave  height  to  length  ratio,  for  sinusoidal,  trochoidal 
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H/L,  WAVE  HEIGHT  TO  LENGTH  RATIO 


Figure  3-3.  Differential  Transmission  of  Orthogonal  Linear 

Polarizations  Averaged  Over  Various  Wave  Spectra. 


and  Stokes  waves.  The  results  show  that  for  values  of  H/L  less  than  the 

■ 2 


Stokes  limit  H/L  of  0.  143  the  differential  polarization  is  less  than  10 
This  means  that  the  effect  at  worst  will  produce  an  error  bias  in  the  Raman 
temperature  measurement  of  about  1°C  or  4 to  8 parts  per  thousand  of  salinity 
concentration.  The  calculation  also  shows  that  the  assumption  of  wave  shape, 
i.  e.  sinusoidal,  trochoidal  or  Stokes,  is  not  a significant  factor  in  the  integral 
of  4B/Bj_  over  the  entire  wave  until  the  Stokes  limit  is  greatly  exceeded. 

If  the  waves  are  small  compared  to  the  laser  spot  size  on  the  water 
surface,  or  if  the  measurement  is  made  by  time  average  of  many  laser 
shots  randomly  distributed  on  large  waves,  the  above  calculations  are  valid. 
For  single  pulse  measurements  with  waves  large  compared  to  the  laser  spot 
size,  the  local  slope  of  the  wave  should  be  monitored  and  a correction  made 
to  the  data.  We  have  made  calculations  that  show  this  to  be  possible. 

3.  1.  1.2  Depth  Resolution.  The  depth  resolution  depends  on  the  laser 
pulse  duration,  the  receiver  bandwidth  and  the  properties  of  the  water 
transmission  medium,  especially  the  refractive  properties  of  the  surface, 
i.  e.  the  air-water  interface.  Off-the-shelf  laser  and  receiver  systems  have 
demonstrated  30  cm  depth  resolution  in  bathymetric  lidar  systems,  v/hich  are 
generically  similar  to  a laser  Raman  range  gated  temperature/ salinity  lidar 
system.  Improvements  in  the  direction  of  shorter  laser  pulse  durations  and 
improved  receiver  bandwidths  and  rise  times  are  possible  to  incorporate  into 
lidar  systems.  The  major  depth  limitation,  at  least  for  an  airborne  system, 
appears  to  be,  however,  the  random  refractive  effects  of  the  air-water 
interface  at  the  sea  surface. 

A computer  routine  similar  to  that  described  in  Section  3.  1.  1.  1 was 
written  to  describe  the  depth  resolution  as  a function  of  sea  state.  The 
listings  of  this  routine  are  exhibited  in  Appendix  C and  is  given  the  name 
SWDCO  (Surface  Wave  Depth  Correction).  Also  exhibited  in  Appendix  C are 
typical  computer  printouts  showing  the  results  of  computer  calculations. 

The  physical  geometrical  optics  basis  for  the  calculation  of  depth 
resolution  as  a function  of  sea  state  is  shown  in  Figure  3-4. 

Figure  3-4  shows  a vertical  laser  beam  entering  a water  surface  that 
is  titled  an  angle  ©i  with  respect  to  the  horizontal.  The  laser  beam  is 
refracted  into  an  angle  with  respect  to  the  normal  © ^ and  attains  an  actual 
depth  that  is  less  than  the  virtual  depth  that  would  have  been  attained  in  the 
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same  transit  time  in  the  absence  of  refraction.  The  fractional  difference 
in  depth  between  the  actual  and  virtual  depth  is  given  by  the  expression 
[ 1 - cos  (©  ^ -©fc)]  . Computer  calculations  were  carried  out,  using  the 
code  listed  in  Appendix  C for  sinusoidal,  trochoidal  and  Stokes  waves  of 
varying  crest  height  to  wavelength  ratio.  The  results  are  shown  in  Figures  3-5 
and  3-6. 

Figure  3-5  presents  plots  of  the  results  of  the  calculations  of  the 
fractional  error  in  depth  measurement  vs.  H/L,  the  peak  to  trough  wave 
height  to  length  ratio.  The  calculations  were  made  for  three  wave  shapes; 
sinusoidal,  trochoidal  and  Stokes.  It  is  apparent  from  the  plots  that  the 
error  in  depth  measurement  is  not  significantly  influenced  by  the  choice  of 
wave  model  for  values  of  H/L  less  than  the  Stokes  limit  of  0.  143.  The  values 
plotted  are  the  fractional  error  in  depth  measurement  as  averaged  over  a 
complete  wave.  Thus,  it  is  seen  that  the  fractional  error  in  depth  measurement 
is  less  than  0.  3%  for  ocean  waves  having  an  H/L  near  the  Stokes  limits  and 
nearer  0.  1%  for  the  more  typical  ocean  waves  with  H/L  from  0.  06  to  0.  08. 

Figure  3-6  presents  three  plots  of  the  results  of  the  calculations  of 
the  fractional  error  in  depth  measurement  vs.  H/L  for  the  4th  order  Stokes 
wave  which  closely  models  real  ocean  waves.  The  plot  labelled  "average" 

presents  the  results  of  averaging  the  fractional  error  in  depth  as  calculated 
at  100  equally  spaced  positions  along  the  wave.  This  is  the  same  curve  as 
plotted  in  Figure  3-5  as  the  Stokes  curve.  The  curve  labelled  "maximum" 
is  a plot  of  the  maximum  fractional  error  in  depth  and  the  curve  labelled 
"90%"  is  a plot  of  the  90th  largest  value  of  fractional  error  in  depth  as 
calculated  over  one  wavelength  of  a wave  for  each  H/L.  It  is  readily  observed 
from  these  plots  that  even  in  the  worst  case  of  the  laser  beam  encountering 
the  maximum  error  condition  (i.  e.  maximum  local  slope)  of  the  wave,  the 
fractional  error  in  depth  measurement  is  less  than  1%  for  waves  having  an 
H/L  near  the  Stokes  limit  of  0.  143.  For  waves  with  smaller  H/L,  more 
typical  of  the  ocean  state,  the  maximum  error  in  depth  measurement  is 
nearer  the  . 1%  to  .2%  range. 

3.  1.  1.3  Lateral  Resolution  The  refractive  effects  of  the  air- 
water  interface  at  the  sea  surface  will  influence  the  lateral  resolution  of  the 


measurement  volume  that  is  attainable  as  a function  of  depth.  Figure  3-4 
shows  the  lateral  error  produced  by  the  surface  refraction.  The  value  of 
the  lateral  error  is  the  product  of  the  depth  and  sin  (0  ^ - © fc). 
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Figure  3-5.  Fractional  Error  In  Depth  Measurement  Due  to 
Refraction  as  Function  of  H/L  For  Stokes, 
Trochoidal  and  Sinusoidal  Waves. 


Fractional  Error  in  Depth  Measurement 


The  relative  importance  of  this  effect,  when  averaged  over  typical 

wave  spectra,  can  be  estimated  using  the  same  computer  program  that  was 

described  in  3.  1.  1.3  above.  As  listed  in  Appendix  C,  on  page  C-7,  the  4th 

order  Stokes  wave  shape,  near  the  Stokes  limit  at  an  H/L  of  0.  1752  exhibits 

an  average  value  of  [ 1 - cos  (0  . -©t)  ] of  5.  08  x 10'3  which  is  the  average 

fractional  depth  error.  Solving  this  expression  for  (0.  - © one  obtains  an 

average  effective  value  for  (©.  -0fc)  of  5.78  degrees.  Substituting  this  value 

into  sin  ( ©.  - ©t)  an  average  value  of  0.  101  is  obtained  for  the  fractional 

lateral  error  for  near  breaking  waves  at  the  Stokes  limit  of  H/'L.  This  means 

that  at  a depth  of  100  meters  the  lateral  position  of  the  measurement  may  be 

in  error  by  as  much  as  + 10  meters  for  Stokes  waves  near  the  limit  of  breaking. 

For  less  violent  waves,  having  smaller  values  of  H/L,  the  error  becomes 

correspondingly  less.  For  example,  again  using  the  data  in  Appendix  C,  on 

page  C-8,  at  a value  of  H/L  = 0.  02  the  average  value  of  l - cos  (0  -©  ) 

„-5  it' 

is  3. 95  x 10  which  gives  an  average  value  of  sin  ( © - © ) of  0.  0089.  This 

means  that  for  waves  with  an  H/L  ratio  of  0.  02  at  a depth  of  100  meters  the 

average  lateral  position  error  due  to  surface  refractive  effects  would  be 

0.89  meters. 

3*1*2  Volume  Transmission  Phenomena  The  effect  of  diffe rential 
attenuation  on  the  application  of  a "two-color"  or  band  shape  analysis  of  the 
liquid  water  Raman  spectrum  has  been  treated  by  previous  investigators. 

The  new  work  on  this  contract  has  extended  the  analysis  to  include  in  addition 
the  effects  of  differential  depolarization  and  the  requirements  for  multiwavelength 
polarization  analysis  for  the  Raman  temperature  and  salinity  measurement. 

3*  1*  2>  1 Differential  Depolarization  - Experimental  The  data 
obtained  by  Duntley'  which  is  shown  in  the  graph  of  Figure  3-7  clearly 
indicates  that  the  decay  of  the  degree  of  polarization  of  polarized  beam  of 
light  as  it  passes  through  natural  waters  is  an  exponential  function  of  the 
distance  traversed.  What  is  remarkable  from  the  Duntley  data  is  that  over 
80%  of  the  original  polarization  remains  even  after  traversing  20  attenuation 
lengths.  Similar  order  of  magnitude  results  for  the  value  of  the  differential 
depolarization  coeffic  vnt  were  obtained  by  other  experimental  investigators. 

In  a recent  private  conversation  with  Duntley  he  expressed  the  opinions 
that  (a)  to  the  extent  investigated,  all  natural  waters  exhibit  the  same  general 
exponential  deploarization  characteristics,  (b)  very  limited  experimentation 
has  been  done  on  the  depolarization  of  transmitted  light  in  a variety  of  natural 
waters,  (c)  the  fluctuations  in  the  depolarization  of  transmitted  light  in  natural 
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Figure  3-7.  Depolarization  During  T ransmission  Through 


Natural  Waters  Obtained  by  S.  Q.  Duntley, 
Reference  11. 


PERCENTAGE  POLARIZATION 


waters  is  virtually  a virgin  territory  for  the  experimentalist  and  (d)  the 
polarization  information  content  of  a beam  of  light  is  preserved  essentially 
until  the  last  photon  has  been  absorbed,  i.  e.  the  attenuation  coefficient  is  very 
much  larger,  the  order  of  100  times  larger,  than  the  depolarization  coefficient. 

The  scaling  of  the  depolarization  during  transmission  through  natural 
waters  and  its  effect  on  a Raman  depolarization  ratio  temperature/  salinity 
measurement  has  been  estimated.  The  data  of  Duntley,  shown  in  Figure  3-7 
justify  the  assumption  of  an  exponential  decay  law  for  depolarization  as 
follows: 

= p j 
dx  S P 

where  1 is  the  polarized  portion  of  a beam  of  light,  x is  the  water  transmission 
distance  and  fl  is  the  depolarization  coefficient.  The  development  of  the 
formalism  in  terms  of  the  observed  depolarization  as  a function  of  the 
depolarization  via  two-way  transmission  through  the  water  as  a noise  source 
in  addition  to  the  water  Raman  scattering  is  given  in  Appendix  D.  The 
result  is: 


P = 1 


1 + p 


P X, 
s 1 


[■('•  X'I 


where  P is  the  observed  ratio  of  orthogonal  polarization  components,  P is 
the  Raman  depolarization  which  is  temperature  sensitive,  Xj  is  the  depth 
from  which  the  scattering  is  being  observed  and  I and  I are  the  amounts 

Uo  Po 

of  upolarized  and  polarized  light  in  the  initial  laser  beam. 

The  results  of  a sample  calculation  for  a typical  case  are  shown  in 
Table  3-1  as  a function  of  KX^  the  depth  in  attenuation  lengths.  The 
parameters  used  as  input  were: 


M 


I.  01  KX, 


= 0.012K  (Reference  2) 

= 0.008K  (Reference  11) 


An  average  value  of  0.01  was  taken  as  the  proportionality  constant 
between  the  diffuse  attenuation  coefficient  K and  the  depolarization  coefficient 
P s,  this  being  an  average  of  the  two  known  literature  values  of  the  constant. 

The  results  in  Table  3-1  show  that  if  the  P effect  were  neglected  and 
not  corrected  for  a temperature  error  AT  of  about  1.6°C  per  diffuse 
attenuation  length  and  a salinity  error  As  of  about  7 parts  per  thousand  per 
diffuse  attenuation  length  could  occur.  The  results  also  show  that  if  the 
proportionality  between  Pg  and  K is  known  to  within  10%  then  after  corrections 
are  made  a residual  uncertainty  of  0.  1°C  in  temperature  and  0.  7 parts  per 
thousand  in  salinity  would  remain  for  each  diffuse  attenuation  length  traversed. 

This  analysis  indicates  that  based  upon  (a)  the  lack  of  experimental 
work  published  or  known  on  the  variability  of  the  depolarization  of  transmitted 
light  in  natural  waters  and  (b)  the  sensitivity  of  the  temperature  and  salinity 
measurement  to  the  P depolarization  coefficient  that  a major  thrust  of  our 
future  research  should  be  to  study  the  extent  of  the  variability  of  the  propor- 
tionality between  P and  K and,  thus,  better  define  the  residual  temperature 
and  salinity  bias  that  can  be  expected. 

The  impact  on  the  experiment  design  and  the  experiment  test  plan  is 
that  a software  capability  to  determine  K through  differential  analysis  of  the 
Raman  signal  as  a function  of  depth  should  be  included.  Independent 
experimental  measurement  of  the  proportionality  constant  between  Pg  and 
K should  be  investigated  for  a variety  of  natural  waters. 

3.  1.2.2  Differential  Depolarization  - Theoretical.  Because  of  the 
importance  of  the  differential  depolarization  effect  described  in  Section  3.  1.2.  1, 
some  effort  was  devoted  to  calculations  to  attempt  to  describe  from  basic  Mie 
scatter  calculations  the  depolarization  in  the  forward  direction  in  natural 
waters.  The  motivation  was  both  to  "explain"  the  previous  experimental  results 
of  relatively  weak,  i.  e.  PgR5  0.01K,  depolarization  in  the  forward  direction 
and  also  to  guide  future  efforts  in  the  directions  of  relfevant  further  calculations 
and/or  experiments.  The  results  of  this  effort  are  summarized  as  follows. 

The  depolarization  of  the  forward- scatte red  component  of  an  under- 
water light  beam  was  investigated.  This  forward-scattered  component 
dominated  the  direct  attenuation  beam  after  a few  scattering  lengths  and  is, 
therefore,  the  significant  component  for  remote  probing  of  sub-surface  ocean 
temperatures.  Computations,  using  the  Mie  formalism  were  made  of  the 
depolarization  within  a five-degree  cone  of  light,  single- scattered  by  spherical 
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particles  which  approximate  in  sizes  and  refractive  indices  those  reported 
for  ocean  water  in  the  literature  . Essentially  no  depolarization  (<<1%) 
was  calculated  for  single  scattering.  This  result  is  in  accord  with  previous 
observations  indicating:  (a)  "negligible"  depolarization  of  the  forward- 
scattered  beam  of  a laser  after  18  scattering  lengths  of  lake  water'  anci 
(b)  20%  depolarization  of  light  after  travelling  20  scattering  lengths  of  lake 
water^^.  Such  polarization  as  occurs  in  natural  waters  may  be  the  result 
of  scattering  from  very  small  particles,  non- sphe r ical  particles  and/or  very 
many  scattering  events.  In  either  case  direct  measurements  of  depolarization 
in  various  types  of  sea-water  would  probably  be  more  useful  than  complex 
and  idealized  calculations. 

Particulates  and  Forward-Scattered  Light  in  Ocean  Water  Brown  and 
Gordon^  ^ have  determined  the  most  likely  size  and  refractive  index 
distributions  of  particulate  matter  in  sea  water  at  a Bahama  Islands  location. 
They  find  the  particulate  matter  can  be  divided  into  three  classes:  (a)  large 
organic  particles  with  a refractive  index  relative  to  water  of  about  1. 03  and 
diameters  ranging  upwards  from  3.  75  urn  with  a sharply  decreasing 
dependence  of  number  density  on  particulate  diameter,  (b)  mid-size  mineral 
particles  with  a relative  refractive  index  of  1.15  and  a size  distribution 
sharply  decreasing  from  1.  25  to  3.  75  n m,  and  (c)  small  organic  particles 
with  relative  refractive  index  of  about  1.01  and  a size  distribution  from  0.65 
to  1.25  fji  m. 

Measurements  and  calculations  by  these  authors  indicate  that  the 
forward- scatte red  beam  was  about  5°  wide  (half  width)  with  a bright  1°  core. 

The  1°  core  was  determined  by  the  large  organic  particles  while  the  beam 
beyond  2°  was  determined  by  the  mid-size  organics.  The  small  organic 
particles  contributed  appreciably  less  than  1%  to  the  forward  scattered  energy. 

The  calculations  for  this  work  were  made  for  particles  with  a diameter 
Dj  of  3.95  m and  a refractive  index  of  1.05  and  for  particles  with  a 
diameter  D ^ of  1.  32  n m and  a refractive  index  m2  of  1.  15.  These  diameters 
were  selected  as  being  representative  since  the  scattering  for  the  large  organics 
and  mid-size  mineral  particles  is  dominated  by  the  smaller  particles  in  each 
group  because  of  the  sharply  decreasing  size  distributions.  Thus  the  calculate  1 
should  give  an  adequate  "feeling"  of  sea-water  depolarization  for  the  exploratory 
purposes  of  this  note. 

Theory  - Large  Particle  Scattering  A theory  has  been  desc ribed  by  Goody^1^ 
for  the  scattering  of  light  by  large  particles  with  (relative)  refractive  index 
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near  unity.  The  theory  should  be  applicable  to  the  large  organic  particles 
discussed  above.  No  depolarization  effects  are  observed  in  this  limit. 
However,  since  the  mid-size  particulate  matter  also  contributes  to  forward- 
scatter  and  estimates  of  small  effects  are  desired,  more  detailed  Mie  theory 
calculations  were  made  as  will  be  discussed  in  the  following. 

Mie  Theory  Computations  The  theory  of  scattering  of  light  by  spherical 
particles  of  arbitrary  size  and  refractive  index  was  developed  by  Mie  and  has 
been  discussed  by  several  authors^14).  We  define  the  "plane  of  observation" 
as  the  plane  containing  the  direction  of  propagation  of  the  incident  wave  and 
the  direction  of  observation.  We  consider  the  intensity  of  light,  (Q  ) of 
plane  polarized  light  whose  incident  electric  vector  is  perpendicular  to  the 
plane  of  observation  and  which  is  scattered  an  angle  ( 0 ) and  plane-polarized 
light,  J ,,  (0)  with  incident  electric  vector  parallel  to  the  plane  of  observation. 
Mie  theory  gives: 
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where  the  A 's  and  B 's  are  complex  functions  of  products  of  Bessel  functions 
n n r 

and  their  derivatives  whose  arguments  depend  upon  particle  properties,  and 

the  P (cos  0 ) are  associated  Legendre  functions.  Both  the  sets  of  A 's  and 
n n 

Bn's  and  the  Legendre  functions  have  been  tabulated  but  by  different  groups^* 
Some  must  be  used  in  computations  using  various  tables  because  of  different 
sign  conventions  for  the  various  quantities. 

For  circularly -polarized  incident  light  the  plane  of  observation  can 
always  be  chosen  so  that  J and  J tl  can  be  considered  as  the  intensities 
along  the  axes  of  the  elliptically  polarized  scattered  light.  The  axis  of  the 
ellipse  will  rotate  with  the  azimuth  of  the  observation  angle. 

Calculations  of  J jl  and  J ,|  were  made  for  © = 1°,  3°  and  5°  for 
the  two  particle  types  (D  = 3.  95  M m,  rr^  = 1. 05,  = 1.  32  u m,  m^  = 1.15) 

discussed  above.  No  depolarization  was  obtained  to  four  significant  figures 
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except  for  the  smaller  particle  size  and  0=5°  for  which  a 0.25%  depolarization 
was  obtained.  (Since  random  depolarization  would  be  expected  to  increase 
something  like  the  square  root  of  the  number  of  scattering  events,  104 
scattering  events  would  be  required  for  a 25%  depolarization.) 

Elaborate  multiple  scattering  calculations  would  be  required  to 
carry  things  much  further.  Also  we  have  neglected  the  scattering  by  the 
smailer  organic  particles  (D  < 1. 25  p m)  which  have  less  than  a 1%  effect 
for  single  scatter. 

In  summary  the  experimental  observations  of  weak  depolarization 
in  the  forward  beam  are  to  be  expected  from  the  character  of  the  Mie  single 
scatter  process  for  the  "large"  particles  which  dominate  sea-water  scatter. 
Direct  experimental  measurement  of  depolarization  in  various  types  of  sea- 
water would  probably  be  more  useful  than  further  complex  and  idealized 
calculations.  We  intend  to  carry  out  such  experiments. 

3 ' 1 ' 1 ’ 3 Differential  Attenuation  and  the  Need  for  MultiwavelenKth 
Polarization  Detection.  The  effect  of  the  well  known!4)  differential  attenuation 
of  light  through  natural  waters  and  its  impact  on  the  Raman  depolarization 
detection  technique  was  also  investigated  as  part  of  the  effort  of  planning  an 
experimental  apparatus  design.  The  main  conclusion  arrived  at  was  that  a 
simultaneous,  multiwavelength  detection  or  wavelength  scanning  detection 
capability  is  required  if  a reasonable  depth  capability  is  to  be  obtained  even 
using  the  polarization  ratio  technique.  (Previous  studies^)  had  suggested 
that  although  differential  attenuation  was  a problem  for  "two  color"  measure- 
ments it  would  not  affect  a polarization  ratio  scheme.  ) 

The  physical  basis  for  this  conclusion  is  that  the  Raman  depolarization 
coefficient  for  liquid  water  varies  across  the  band,  ranging  for  circular 
polarization  from  values  of  0.  2 to  0.  8.  If  the  entire  Raman  band  is  detected 

with  a single  detector  an  average  P will  result. 

R 

The  average  P R is  effectively, 

MX)  P R ( X ) dX 
R Jl  (X)  dx 

where  the  PR  ( X ) is  weighted  by  the  intensity  I ' ( X ). 
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Differential  attenuation,  as  for  example  would  be  caused  by  traversing 
several  attenuation  lengths  of  water,  would  cause  a distortion  of  the  Raman 
band  shape  i.  e.  the  I ( X)  factor,  and  thus  introduce  an  error  in  the  effective 
P ^ measured.  A solution  is  to  use  multiple  band  segments  of  a bandwidth 
such  that  the  percentage  distortion  is  acceptable  in  each  band  segment.  The 
polarization  ratios  are  thus  obtained  on  a band  segment  by  band  segment  basis 
and  then  combined  with  suitable  weighting  functions  to  yield  the  temperature 
and  salinity  information. 

The  question  addressed  in  the  present  new  work  was  the  relationship 
between  the  bandwidth  of  the  individual  Raman  band  segments,  the  accuracy 
required  and  the  differential  attenuation.  The  practical  result  was  that  a 
required  spectrometer  resolving  power  could  be  specified  for  the  experimental 
apparatus  design. 

Calculations  were  carried  out  using  the  diffuse  attenuation  data  of 
(4) 

Jerlov  for  natural  ocean  water  which  can  be  fit  by  a parabolic  function 

2 

such  as  K = /?  ( X , - X ) +7  where  0 and  7 are  constants  which  are  functions 

^ -6-1-2 
of  water  quality  and  type.  Using  a typical  value  of  0 of  4 x 10  m nm  ■ 

results  were  obtained  showing  the  residual  error  remaining  as  the  Raman 

band  was  piecewise  detected  by  means  of  multiwavelength  analysis. 

Based  on  this  analysis,  the  fractional  change  in  signal  between  two 

wavelengths  X ^ and  X^  caused  by  differential  transmission  compared  to  a 

uniformly  transparent  medium  is  given  by  the  following  equations. 


1 ACTUAL 
IUNIFORM 
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Uniform 
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The  plots  in  Figures  3-8  and  3-9,  based  on  the  above  equations  show 
the  signal  error  introduced  as  a function  of  the  bandwidth  with  depth  as  a 
parameter.  Figure  3-8  is  for  Xj  - XQ  = 0,  meaning  symmetric  alignment  of 
the  Raman  band  with  the  minimum  of  the  diffuse  attenuation  curve  and 
Figure  3-9  is  for  Xj  - XQ  = 10  nm,  i.  e.  10  nm  displacement  of  the  Raman 
band  from  the  minimum  of  the  K curves. 


25 


! 


0.1  0.2  0.4  0.8  2.0  4.0  8.0 


X2  - A.,  (NANOMETERS) 

Figure  3-9.  Differential  Attenuation  vs . Wavelength 
Separation  With  Depth  as  a Parameter 
(100  £ From  K minimum). 


Based  on  these  calculations  it  is  clear  that  an  experimental  wavelength 
scanning  capability  with  a wavelength  resolution  the  order  of  a few  tenths  of 
nanometers  is  necessary.  It  is  also  obvious  from  these  calculations  that  there 
is  a requirement  for  a tunable  laser  to  position  the  Raman  water  band  at  or 
near  the  minimum  of  the  diffuse  attenuation  (K)  vs.  wavelength  curve.  Both  of 
these  considerations  were  factored  into  the  experimental  apparatus  design  plan. 

3.1.3  Non-Raman  Photon  Sources.  Laser  induced  fluorescence 
backgrounds  as  well  as  ambient  light  backgrounds  are  common  sources  of 
interference  in  remote  Raman  measurement  systems.  Such  effects  can  be 
reduced  by  using  the  following  standard  techniques. 

The  detection  system  should  be  limited  to  only  those  wavelengths  that 
are  within  the  Raman  band  and  synchronous  time  gating  techniques  should  be 
employed.  The  use  of  background  subtraction  by  sampling  the  ambient  back- 
ground during  the  laser  interpulse  period  should  also  be  included  in  the 
experiment.  In  addition,  subtraction  of  laser  induced  fluorescence  can  be 
accomplished  by  measurement  of  the  laser  induced  spectral  background  level 
at  wavelengths  adjacent  to  the  Raman  water  band.  This  capability  should 
also  be  provided  in  the  experiment  apparatus  design. 

The  following  discusses  the  salient  features  of  ambient  light  and  laser 
induced  fluorescence-like  phenomena  when  viewed  as  non-Raman  photon  noi.se 
sources . 

3.  1.3.1  Ambient  Light  Background.  Standard  references^  ^ on 
remote  laser  Raman  applications  show  that  the  limting  signal  to  noise  ratio 
that  can  be  achieved  in  a remote  Raman  system  is  the  "shot  noise"  in  the  laser 
Raman  signal  itself.  Such  a Raman  system  is  assumed  in  the  planned 
expe  riments. 

The  contribution  of  ambient  background  illumination  to  the  signal  to 
(18) 

noise  ratio  can  be  shown'  ' to  increase  the  time  required  to  make  a Raman 
measurement  with  a fixed  signal  to  noise  ratio  as  follows: 


(— ) — L 

V N / P A1 


NEP 


where  T isthe  time  required  to  make  a measurement,  S, N is  the  signal  to  noise 
ratio,  ^AyQ  is  the  average  laser  power,  Pp^,  is  the  peak  laser  power  and  Pj^p 
is  the  noise  equivalent  power  of  the  ambient  background  illumination. 
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Thus  the  important  parameter  is  the  ratio  of  the  background  illumination  or 

the  "noise- equivalent-power"  t'o  the  laser  peak  power. 

The  upward  irradiance  near  the  sea  surface  is  given  by  Jerlov^  as 
2 

being  7 mW /cm  - n m at  a wavelength  of  500  nanometers.  Using  this  number 

for  irradiance  and  neglecting  the  effects  of  sun  glints  and  specular  reflection 

a to  rati°  slightly  greater  than  unity  was  computed  for  an  airborne 

platform  operating  at  an  altitude  of  300  meters  and  using  the  Raman  system 

parameters  assumed  in  the  calculation  of  Appendix  A. 

In  a boat  borne  experiment,  because  the  Raman  signal  level  using  the 

same  laser  source  will  be  considerably  higher,  the  ratio  is 

N ha  x j Jr  x\. 

expected  to  be  small  compared  to  unity.  However,  because  better  than  1 % 
measurements  are  desired  a capability  should  be  designed  into  the  experiment 
to  sample  the  background  ambient  light  level  in  the  time  interval  between 
laser  pulses.  This,  background  measurement  should  independently  sample 
each  of  the  two  polarization  channels  so  that  background  subtraction  may  be 
applied  independently  to  each  polarization  signal  before  a polarization  ratio 
is  computed.  This  is  important  because  ambient  light  from  sea  water  is 
known  to  be  partially  polarized. 

3.  1.3.  2 Laser  Induced  Fluorescence.  The  presence  of  laser  induced 

(18) 

fluorescence-like  background  noise  has  been  shown  ; to  add  another  term  to 
the  equation  which  expresses  the  time  required  to  make  a Raman  measurement 
with  a fixed  signal  to  noise  ratio  as  follows: 


where  the  additional  term  F is  defined  as  the  laser  induced  background 
fluorescence  level  to  Raman  signal  ratio. 

New  experiments  in  natural  waters  should  be  designed  to  directly 
measure  the  laser  induced  fluorescence  to  Raman  ratio  F so  that  estimates 
of  measurement  times  with  constant  signal  to  noise  ratio  can  be  made  as  a 
function  of  laser  excitation  wavelength.  This  capability  is  recommended  as 
part  of  the  experiment  apparatus  design. 

It  should  be  noted  that  preliminary  field  experiments  that  were 
conducted  in  natural  coastal  waters  by  CGC  personnel  and  which  are  described 
in  Section  3.2.3  of  this  report  showed  a fluorescence  to  Raman  ratio  the  order 
of  0.03  or  less  when  using  a 337  nanometer  laser  source.  Such  data  should  be 
extended  to  cover  the  entire  visible  spectrum  by  using  a tunable  laser  source. 
Such  spectral  coverage  should  be  obtained  from  a variety  of  coastal  water  types. 
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3.  2 Detailed  Experimental  Apparatus  Design 

This  section  describes  the  detailed  experimental  apparatus  design 
that  has  been  developed  as  a means  of  carrying  out  experimental  feasibility 
studies  on  the  use  of  laser  Raman  spectroscopy  for  the  remote  measurement 
of  subsurface  temperature  and  salinity.  Section  3.2.  1 describes  a laser 
Raman  spectrometer  that  would  be  appropriate  for  conducting  the  studies; 

Section  3.2.2  describes  the  CGC  research  vessel  and  the  equipment  that  is 
currently  available;  and  Section  3.  2.  3 describes  preliminary  experimental 
data  that  has  been  obtained  from  natural  waters  using  the  CGC  research  vessel. 

3.2.1  Laser  Raman  Spectrometer  for  Field  Utilization 

The  functional  block  diagram  of  Figure  3-10  illustrates  an 
experimental  laser  Raman  measurement  system  that  will  enable  experiments 
to  be  conducted  so  that  the  important  questions  concerning  interferences  and 
noise  sources  can  be  answered. 

t 

The  Government  owned  laser  that  is  currently  available  at  CGC  for 
use  on  this  program  is  a pulsed  nitrogen  laser  capable  of  producing  100  kilowatt 
peak  power  pulses  of  10  nanoseconds  duration  at  a pulse  repetition  rate  of 
500  pulses  per  second.  The  wavelength  of  the  laser  is  337  nanometers. 

The  output  from>the  nitrogen  laser  is  used  to  pump  a dye  laser 
which  can  be  tuned  over  thje  entire  range  of  interest  in  the  blue-green  visible 
spectrum,  i.  e.  450  to  500  nafiometers.  The  output  from  the  dye  laser  is 
conditioned  by  sequential  passage  through  an  interference  filter  to  block  stray 
light  at  the  Raman  frequencies  and  through  a linear  polarizer  and  a quarter 
wave  plate  to  produce  a pure  circularly  polarized  laser  beam.  This  beam  is 
directed  at  the  water  volume  to  be  measured. 

The  Raman  backscattered  light  is  collected  by  a dual  polarization 
receiver  in  which  right  and  left  circularly  polarized  light  is  separately  collected 
and  wavelength  analyzed.  A dual  1-meter  spectrometer  modified  to  operate 
in  a parallel  mode  rather  than  in  a series  mode  is  shown  as  the  wavelength 
analyzer.  This  would  provide  the  few  tenths  of  a nanometer  resolution  that  the 
differential  attenuation  and  polarization  ratio  detection  study  showed  was 
necessary  and  which  was  described  in  Section  3.  1.2.3. 

As  shown  in  the  block  diagram  the  dual  output  from  the  spectrometer 
is  detected  by  a pair  of  photomultipliers,  PMi  and  PM2  in  Figure  3-10,  the 
output  of  which  is  sent  to  a pair  of  coincidence  circuits.  The  other  input  to  the 
coincidence  circuits  is  a gate,  suitably  delayed  relative  to  the  laser  firing,  which 
determines  the  depth  from  which  the  Raman  scattering  is  being  detected.  This 
photon  counting  mode  of  operation  in  which  either  a "0"  or  a "1"  is  detected 
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Figure  3-10  Functional  Block  Diagram  of  Dual  Polarization  Raman  Spectrometer 


for  each  laser  pulse  has  been  shown  to  be  highly  accurate  and  linear  method 
of  conducting  exploratory  laser  Raman  experiments. 

The  question  of  maintaining  the  "balance"  of  two  photomultipliers 
is  always  important.  In  a "0"  or  "1"  mode  of  operation  the  constancy  of  the 
photoefficiency  of  the  photocathode  surface  is  the  important  factor,  not  the 
gain  of  the  dynode  chain.  In  addition  the  two  polarization  channels  can  be 
optically  alternated  relative  to  the  two  photomultipliers  to  enable  cross 
calibrations  on  a regular  basis. 

The  direct  output  from  the  dual  polarization  coincidence  circuits 
is  displayed  to  enable  maximum  experimentalist/operator  interaction  in 
real  time.  A minicomputer  processor  (CPU)  is  used  to  provide  both  real  time 
control  functions  as  well  as  data  acquisition  and  recording.  The  bLock  diagram  of 
Figure  3-10  illustrates  the  experimental  configuration  showing  the  computer 
with  the  spectrometer  control  function  and  the  coincidence  circuit  data  flow 
passing  through  the  processor  input/output  interface  (I/O). 

The  computer  oriented  hardware  configuration  shown  in  the  block 
diagram  of  Figure  3-10  will  accommodate  several  software  approaches  to 
maximize  the  information  content  of  the  data  that  is  acquired.  Examples  of 
two  such  operational  procedures  are: 

a)  By  sweeping  many  times  through  the  spectrum  of  interest 
with  relatively  short  data  taking  periods  at  each  wavelength, 
the  errors  introduced  by  drift  of  various  system  parameters 
may  be  minimized.  Such  a scanning  procedure  is  readily 
permitted  by  the  computer  controlled  spectrometer. 

b)  Scan  the  spectrum  of  interest  in  a manner  which  yields  a 
constant  signal-to-noise  ratio  for  the  data.  This  is 
accomplished  by  scanning  the  spectrum  at  a rate  proportional 
to  the  strength  of  the  signal  being  detected. 

3.2.2  Research  Vessel  Facility* 

A basic  premise  of  this  program  is  to  substantially  advance  the 
techniques  of  water  Raman  spectroscopy  in  situ.  As  such,all  experiments/ 
tasks  (paragraph  3.  3.  1)  are  to  be  implemented  as  a function  of  variations  in 
natural  water  qualities.  The  options  include;  further  laboratory  tank  tests  or 
operating  from  a pier  or  boat. 

It  is  difficult  if  not  impossible,  to  adequately  simulate  natural 
coastal  water  properties  and  gradients  in  tanks.  Important  optical  properties 
of  the  water,  i.  e.  depolarization  properties,  fluorescence  characteristics  and 

* 

The  work  described  in  Section  3.  2.  2 was  supported  in  part  by  NASA  Wallops 
Flight  Center. 
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its  turbidity  depend  to  a large  extent  upon  living  biological  organisms,  many 
of  which  cannot  sustain  life  in  tanks.  It  is  important  that  the  activities  have 
access  to  freshwater  intrusion,  natural  sediment  loading  and 'acquatic' life 
and  temperature  gradients  to  directly  address  coastal  applications. 

A survey  of  pier  and  bridge  locations  from  nautical  charts  as 
well  as  practical  considerations  such  as  the  availability  of  power  and 
coordinating  traffic  patterns  rapidly  leads  to  the  conclusion  that  a completely 
self  contained  floating  platform  is  the  more  cost  effective  way  to  access  the 
necessary  operational  environments. 

The  requirements  and  characteristics  of  a small,  economical  and  easy 
to  operate  oceanographic  research  vessel  have  been  exhaustively  investigated 
with  internal  support.  These  files  are  available  in  our  Wakefield  office. 

A 35'  twin  screw,  nine  ton  displacement  research  vessel,  Makai, 

9jt 

Figure  3-11  has  been  fitted  out  and  currently  hosts  the  operational  apparatus 
configuration  shown  as  Figure  3-12  which  was  used  to  obtain  the  first  in  situ 

data  discussed  under  Section  3.  2.  3. 

The  current  laser  Raman  system  is  a bistatic  configuration.  The 
laser  source  is  a pulsed  nitrogen  laser  operating  at  a wavelength  of  337.  1 
nanometers  and  producing  pulses  of  100  kilowatts  peak  power  with  an  effective 
pulse  duration  of  10  nanoseconds  at  a pulse  repetition  rate  of  500  Hertz. 

As  shown  in  t>he  block  diagram,  a bistatic  arrangement  of  the 
transmitter  and  receiver  optics  function  to  define  a volume  of  subsurface 
water  from  which  the  laser  Raman  scattering  is  detected.  Operation  in  such 
a bistatic  mode  is  adequate  for  experimental  studies  in  the  first  few  tens  of 
meters  of  depth.  Practical  application  of  the  technique  should  utilize  electronic 
range  gating  techniques,  especially  for  deeper  penetration  applications.  The 
photograph  of  Figure  3-13  shows  the  equipment  in  the  laboratory  undergoing 
calibration  prior  to  installation  aboard  the  Makai. 

The  photons  collected  by  the  receiver  optics  are  passed  through  a 
double  1/4  meter  focal  length  spectrometer  having  a 0.  5 nanometer  spectral 
resolution.  The  function  of  the  spectrometer  is  to  produce  a spectral  scan 
of  the  liquid  water  Raman  O-H  stretching  band  with  sufficient  resolution 
so  that  the  monomer  and  polymer  components  can  be  resolved  and  a temperature 
determined. 


■?  ‘ 

Much  of  the  Government  owned  equipment  is  accountable  under  NASA 
Contract  NAS6-2751 
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Figure  3-12 


Block  Diagram  of  Laser  Raman  System  Currently 
Operational  on  Research  Vessel  Makai. 
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The  Raman  photon  output  from  the  spectrometer  is  detected  by 
a RCA  1P28  photomultiplier,  amplified  by  a HP  461  A amplifier  and  conditioned 
by  a HP  pulse  electronics.  The  conditioned  Raman  photon  signal  and  a gate 
pulse  suitably  synchronized  and  delayed  with  respect  to  the  laser  firing  were 
combined  in  a coincidence  circuit  the  output  from  which  is  recorded  by 
counter  #1.  The  number  of  laser  pulses  occurring  is  recorded  in  counter  #2. 

The  ratio  of  the  value  recorded  by  counter  #1  to  counter  #2  is  the  Raman 
signal  obtained  over  a given  time  interval. 

The  photograph  of  Figure  3-14  provides  a close-up  view  of  the  optical 
transceiver  cart  configuration  selected  for  implementation.  A custom 
constructed  mobile  cart  houses  the  laser  transmitter  and  spectrometer 
receiver  in  an  over  and  under  configuration.  Rigid  extensions  of  the  shelves 
are  used  to  support  the  transmitter /receiver  optics  and  extend  the  "viewing 
points"  45cm  beyond  the  transom  of  the  vessel  where  high  quality  folding  mirrors 
are  used  to  complete  the  path. 

The  compact  120cm(L)  x 75cm(W)  x 97cm(H)  cart  is  equipped  with 
wheels  which  facilitate  service  access  and  withdrawal  into  the  cockpit  for 
protected  storage  when  not  in  use. 

The  laser  employed  (top  shelf)  is  a GFP  Avco  C950  nitrogen  unit 
with  a CGC  modified  cooling  system  in  order  to  lighten  payload.  The 
spectrometer  (bottom  shelf)  is  a Spex  model  1672  Doublemate  and  1200  g/mm 
grating  set  blazed  at  300  nm. 

3.2.3  Preliminary  Experimental  Raman  Data  Obtained  from 
Natural  Waters** 

The  Raman  lidar  system  that  is  shown  in  Figure  3-12  and  is  now 
currently  on  the  research  vessel  Makai  has  been  operated  to  obtain  preliminary 
laser  Raman  spectra  from  in  situ  natural  coastal  waters.  The  data  was 
analyzed  to  obtain  a temperature  measurement.  It  is  important  to  note  that 
this  is  the  first  such  data  known  to  have  been  taken  in  natural  coastal  waters. 

The  Raman  spectra  were  obtained  from  the  Annisquam  River  in 
Massachusetts  near  the  berth  of  the  research  vessel.  Spectra  were  obtained 
over  nearly  a full  tidal  cycle  with  no  apparent  change  in  the  character  of  the 


Typical  spectral  data  that  were  obtained  in  the  Annisquam  River  are 

shown  in  Figure  3-15.  This  i’s  a plot  of  the  number  of  Raman  photoeiectrons 

collected  per  0.  5 nanometer  spectral  interval  per  10,  000  laser  pulses  as 

a function  of  wavelength.  Shown  are  two  sets  of  data,  the  dashed  curve 

obtained  as  an  instrument  calibration  in  the  laboratory  with  a water  sample 

at  22° C and  the  solid  curve  obtained  in  the  Annisquam  River  water  at  6°C. 

The  temperature  is  obtained  by  normalizing  the  peaks  and  observing  the 

change  in  height  of  the  band  contour.  A temperature  scale  has  been 

calculated  based  on  laboratory  data  and  placed  on  Figure  3-15  so  that 

22 °C  corresponds  to  the  laboratory  calibration  data.  The  Annisquam 

River  temperature  of  4°C  inferred  from  the  spectra  is  consistent  with 

the  measured  value  of  6°C  to  within  the  statistical  precision  of  the 

measurement  which  was  + 2 C.  Large  improvements  in  accuracy  will 

be  observed  by  increasing  the  photoelectron  sample  size  and  also  by 

(19) 

using  standard  least  square  computer  curve  fitting  to  enable  analysis 
of  the  entire  Raman  band. 

A short  research  note  has  been  written  and  published  in 
Geophysics  Research  Letters  to  describe  these  initial  results.  A preprint 
of  the  note  is  included  as  Appendix  E of  this  report. 
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Figure  3-15  Raman  Spectral  Data  for  Liquid  Water. 
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Experimental  Field  Test  Plan 

A series  of  field  experiments  in  a variety  of  natural  coastal  water 
conditions  are  required  to  further  study  and  quantitatively  determine  the 
extent  of  natural  and  equipment  related  interferences  and  noise  sources 
(Paragraph  3.  1)  upon  the  laser  Raman  technique  to  measure  subsurface 
temperature  and  salinity. 

Abroad  plan  for  Raman  subsurface  temperature /salinity  research 
has  been  considered  to  carry  the  program  from  its  present  formative  stage  to 
a more  developed  and  predictable  documented  technique.  The  Schedule  and 
Plan  for  Raman  Tempe  rature /Salinity  Research  Figure  3-16  illustrates  a 
modestly  supported  (1-  1 1/2  man  years/yr)  program  which  accommodates 
10-12  month  period  of  primary  experimentation  and  apparatus  improvement. 

During  the  second  phase  it  is  anticipated  that  the  refined  and 
developed  experimental  apparatus  will  be  capable  of  processing  large 
volumes  of  data  in  a real  time  profiling  mode  through  substantial  temperature 
and  salinity  gradients.  It  is  desirable  to  validate  and  further  develop  the 
techniques  in  the  presence  of  more  accurate  and  sophisticated  sea  truth 
instrumentation,  preferably  as  an  adjunct  to  an  already  established  field  research 
or  monitoring  program.  Several  candidate  locations  capable  of  supporting 
such  a coordinated  effort  have  been  considered  and  they  are  discussed 
under  Paragraph  3.  3.  2.  Sea  Truth  Requirements. 

The  third  phase  which  will  provide  a broad  data  base  in  deep  ocean 
as  well  as  coastal  waters  might  very  well  take  the  form  of  installing  the 
apparatus  aboard  one  or  more  of  the  manned  vessels  operated  under  the  support 
of  the  Office  of  the  Oceanographer  for  research  and  data  base  gathering  purposes. 
One  such  vessel  which  in  our  view  is  particularly  well  suited  to  the  task  is  the 
twin  hulled  U.  S.  N.  S.  HAYES. 

Of  the  specific  initial  experimental  tasks  listed  below  in  Paragraph 
3.  3.  1 some  must  be  implemented  in  series.  All  experiments  are  as  a function 
of  natural  water  quality  (i.  e.  diffuse  attenuation  coefficient  (K),  fluorescence 
to  Raman  ratio  (F),  salinity,  sediment  loading,  aquatic  life,  etc.  ) and  depth. 
Laboratory  experiments  have  been  performed  by  Walrafen,  Schwiesow,  ^2  1 ^ 
Chang  and  Young,  ^ ^ and  Slusher  and  Derr.  However,  it  can  be  readily 
seen  that  it  is  diffucult  if  not  impossible,  to  simulate  natural  conditions  in 
laboratory  tanks. 
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SCHEDULE  FOR  RAMAN  TEMPERATURE/SALINITY  RESEARCH 


r 


The  most  cost  effective  means  to  access  natural  water  quality 
conditions  is  to  place  an  apparatus  aboard  a boat  capable  of  relocating. 

The  apparatus  should  be  designed  so  as  to  be  transportable  by  two  persons 
by  truck  to  boats  operating  out  of  different  geographical  locations.  A 
maximum  gross  apparatus  weight  under  900  kg  with  power  consumption 
under  6KW  will  provide  for  accommodation  aboard  average  vessels  9 meters 
long  and  greater.  Afunctional  block  diagram  of  an  existing  Government 
owned  apparatus  available  to  support  this  field  test  plan  appears  as  Figure  3-12. 
The  preferred  configuration  which  can  be  readily  assembled  from  additional 
existing  Government  inventory  and  minor  purchases  is  shown  as  Figure  3-10. 

The  balance  of  this  section  discusses  experimental  objectives  and 
methods,  sea  truth  requirements,  general  test  site  considerations  and  a well 
researched  and  documented  initial  recommended  testing  region  with  supporting 
nautical  charts. 

3.3.1  Experimental  Objectives  and  Method 

3.  3.  1.  1 Transmission  (i,  e,  K)  Derive  and  catalog  a 
class  of  attenuation  coefficients  as  a function  of  various  natural  coastal  water 
qualities  and  as  a function  of  depth.  Operate  a tunable  dye  laser  through  the 
wavelength  regions  400  to  500  nanometers.  Observe  the  Raman  shifted  spectral 
region  between  420  and  580  nanometers.  Perform  slope  analysis.  Correlate 
with  alpha  meter  readings.  Evaluate  data  to  specify  optimal  1/4  wave  plate 
for  advanced  experiments. 

3.  3.  1.2  Fluorescence  Raman  Ratio  (F)  Survey  the  Raman 
to  Fluorescence  ratio  as  a function  of  various  natural  coastal  water  quality. 
Operate  tunable  dye  laser  as  in  3.  3.  1.  1.  Exercise  as  a function  of  depth  in 
several  varying  water  quality  locations  and  establish  an  historical  data  base. 

3.3.  1.3  Differential  Polarization  Analysis  Operate  the 
dye  laser  in  selected  spectral  regions  based  upon  the  results  of  3.  3.  1.  1 and 

3.3.  1.2  above.  Obtain  spectrally  resolved  polarization  ratio  data.  Evaluate 
data  to  verify  required  spectral  resolution.  Specify  multichannel  wavelength 
detection  requirements. 

3.3.  1.4  Proportionality  Between  p and  K.  Analyze  the 

s 

data  obtained  in  3.3.  1. 3 in  accordance  with  the  scaling  laws  described  in 

Appendix  D and  discussed  in  3.  1.2.  1.  Determine  the  variability  in  the 

proportionality  between  P and  K as  a function  of  the  type  of  water,  water 

s 

quality  and  the  absolute  levels  of  p and  K. 

s 
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3.3.  1.5  Temperature  vs.  Salinity  Separability  The 
spectrally  resolved  polarization  ratio  data  obtained  in  3.3.  1.  3 will  also 
be  analyzed  to  show  the  extent  to  which  temperature  and  salinity  may  be 
independently  measured. 

3.3.2  Sea  Truth  Requirements 

Reliable  and  accurate  sea  truth  and  meteorological  data  is  required. 
The  options  to  obtain  sea  truth  range  from  totally  R/V  self  contained 
instrumentation  to  heavy  dependence  upon  coordinating  experiments  ("piggy- 
backing") with  existing  hydrographic  activities.  Figure  3 - 1 6 shows  the 
desirable  schedule  and  increasing  experiment  diversification  plan.  During 
the  initial  Gulf  of  Maine  field  effort  it  is  felt  that  very  modest  sea  truth 
instrumentation  obtainable  from  Government  inventory  or  local  purchase 
will  be  adequate. 

As  discussed  earlier  in  this  section,  the  second  phase  will  be 
oriented  towards  a greater  dynamic  range  of  performance  and  data  base 
gathering,  suggesting  more  sophisticated  automatic  data  logging  equipment. 
Towards  this  end  and  as  part  of  the  planning  effort  we  have  been  in  contact 
with  cognizant  personnel  at  NASA  Wallops  Flight  Center,  NOAA  Remote  Sensing 
Laboratory,  Miami,  and  Department  of  Natural  Resources  - State  of  Maryland 
(DNR).  All  have  expressed  a strong  interest  in  promoting  the  state  of  the  art 
by  making  available  space  on  their  own  vessels  of  opportunity  at  little  or  no 
cost  to  this  program. 

NOAA's  east  coast  activities  are  supported  by  three  research 
vessels;  a 83m  state-of-the-art  instrumentation  equipped  vessel  capable 
of  open  ocean  ground  truth  measurements,  a 45m  instrumented  vessel 
operating  in  and  about  the  port  of  New  York  and  a smaller  high  speed  vessel. 

Sea  truth  observing  techniques  include:  continuous  subsurface 
temperature  devices,  continuous  salinity  monitors,  optical  spectrometers, 

Sechi  disc,  expendable  bathy  thermographs  (XBT),  drifting  and  moored 
buoys,  etc. 

The  DNR-Maryland  facility  under  Dr.  Myron  Miller  is  responsible 
for  the  State  of  Maryland  Power  Plant  siting  program.  A regular  function  of 
this  office  is  to  monitor/measure  thermal  plume  discharges  in  three  dimensions 
from  power  plants.  An  offer  has  been  extended  to  place  the  apparatus  aboard  a 
DNR  vessel  where  additional  sea  truth  will  be  obtained  from  their  highly 
resolved  thermistor  chain  which  is  towed  through  the  plume  in  grid  pattern 
and  well  documented  during  three  seasons/yr. 
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Additional  candidate  sites  and  coordination  planning  for  Phase  3 
is  better  left  for  a later  point  in  this  program. 

3.  3.  2.  1 Sea  Truth  Measurement  Parameters  At  a 
minimum  records  for  all  phases  should  be  kept  to  include: 

a)  date,  time,  investigator's  name 

b)  location  of  experiment  (latitude  and  longitude) 

c)  sea  state  and  surface  conditions  (dialog  and  photos) 

* d)  atmospheric  temperature,  humidity,  wind 

e)  sea  surface  and  subsurface  temperatures 

f)  salinity  (Initially  surface  water  samples  will  be 
drawn  at  each  location  and  portions  of  water  will  be 
sent  to  the  laboratory  for  analysis.  The  second 
portion  will  be  tagged  by  experiment  number  and 
preserved.  Salinity  gradients  may  be  tested  through 
the  use  of  classical  Nansen  jars.  ) 

In  addition  to  the  above  the  computer  will  automatically 
calculate  and  tabulate  K,  the  diffuse  attenuation  coefficient,  and  the 
background  fluorescence  level  in  the  wings  of  the  Raman  spectra. 
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Test  Site  Requirements 


Proper  selection  of  experimental  test  sites  is  a vital  ingredient 
toward  the  development  and  validation  of  new  remote  measurement  and 
refinement  techniques.  Selected  sites  must  offer  a wide  variety  of  water 
quality  conditions  including  a range  of  attenuation  coefficients, representative 
cross  section  of  coastal  and  estuary  impurities  (fluorescence  sources), 
diverse  micro  organism  t/pe  count.and  a range  of  salinity  conditions.  While 
variations  in  sea  floor  types  and  bottom  vegetation  is  considered  to  have  a 
minimal  effect  on  a gated  Lidar  system  variation  in  texture,  grain  size, 
density  and  depth  in  the  presence  of  tidal  current  or  wave  turbulance  could 
produce  significant  perturbations  in  turbidity  patterns  on  a time  scale  which 
is  short  relative  to  an  experimental  period.  The  natural  vicissitudes  of  the 
ocean  and  veather  conditions  can  be  counted  upon  to  offer  wide  variations 
in  sea  surface  conditions. 


Within  reasonable  distances  of  the  R/V  cruising  range  it  is 
desirable  to  access  zones  of  relatively  deep  clear,  shallow  clear,  deep- 
turbid,  and  shallow  turbid  water  which  are  expected  to  vary  synoptically 
and  in  a reasonably  predictable  pattern. 


3.3.1  Water  Qualitx 


Of  the  preceding  site  requirements  water  quality  is  the  factor 


of  initial  major  experimental  interest. 


Turbidity  in  sea  water  is  due  primarily  to  the  presence  of 
tiny  particulate  matter  or  plankton  entrained  in  the  water.  The  latter 
source  can  be  maximized  or  minimized  by  selecting  winter  or  summer 
frames  of  operation.  Increases  in  turbidity  are  generally  a result  of 
the  intrusion  of  suspended  particulate  matter  into  a watermass  from  an 
extraneous  source  (such  as  runoff  from  land  following  percipitation)  or 
the  resuspension  of  fine  bottom  materials  as  a result  of  the  action  of  tidal 
currents  or  wave  turbulence.  In  any  case,  these  increases  in  particulate 
matter  are  usually  rapid,  with  pronounced  changes  in  turbidity  being  measured 
in  hours.  Decreases  in  turbidity,  on  the  other  hand,  are  usually  slow  to  occur, 
often  requiring  days  or  even  weeks  for  a water  mass  to  return  to  a background 
level.  A downward  flux  in  turbidity  in  an  area  of  water  results  from  the 
settling  of  the  entrained  particulate  matter  or  from  the  replacement  of  the 
watermass  by  clearer  water.  The  former  process  is  extremely  slow  and 
requires  a period  of  relatively  low  turbulence.  The  latter  process  requires  a 
source  watermass  and  a constant  current  or  drift  from  the  clear  watermass 


into  the  turbid  watermass. 


It  has  been  well  .documented  that  salinity  concentrations  in  open 


ocean  vary  slowly  over  hundreds  of  square  miles.  In  addition  open  ocean 
surface  layers  are  usually  well  mixed  often  leading  to  the  assumption  of 
salinity  as  a constant. 

Unlike  ocean  waters,  salinity  in  coastal  waters  vary  dramatically 
responding  to  the  intrusion  of  fresh  water  from  rivers,  heavy  rain  runoffs, 
snow  and  ice  melt  runoff  along  the  coast.  In  the  summer  during  long  periods 
of  zero  rainfall  and  above  normal  temperatures  significant  increases  in 
salinity  may  be  expected. 

Performing  experiments  in  spring  through  late  fall  in  coastal 

waters,  latitudes  40°  and  greater  offer  a range  of  salinity  concentrations 

(22) 

from  ~ 9 ppt  to  in  excess  of  30  ppt,  while  navigable  access  to  the  head 
waters  of  a major  drainage  basin  will  offer  opportunities  to  data  sample 
nearly  fresh  waters. 

Initial  Test  Site 

Cape  Ann  Marina  - Experimental  apparatus  shakedown  and 
initial  calibrations  begin  in  the  laboratory  and  are  refined  in  the  field 
establishing  a feasibility  data  base  of  initial  equipment  performance  in 
natural  waters.  The  waters  surrounding  Cape  Ann  from  41°  10'  to 
42°  50'  latitude  and  71°  02'  to  71°00',  Figure  3-17  provides  suitable  water 
quality  and  sea  conditions  with  convenient  access  for  this  early  task 
including  the  clear  water  masses  in  and  around  Annisquam,  fresh  water 
intruding  from  the  Merrimack  River  and  the  turbid  waters  of  Boston 
Harbor.  Access  to  Power  Plant  thermal  plumes  are  available  in  Salem  and 
Boston  Harbors. 

The  oceanographic  characteristics  of  this  region  are  well 
documented  and  have  been  the  subject  of  numerous  studies  supported  by 
UNH  Sea  Grant  - Raytheon  Program,  NOAA  Environmental  Research 
Laboratory  and  Massachusetts  Department  of  Natural  Resources.  Much 
of  the  characterizations  presented  here  were  documented  in  the  Proceedings 
of  Conference  on  Toxic  Dinoflag ellate  Blooms  hosted  by  Massachusetts 
Science  and  Technology  Foundation  and  Massachusetts  Institute  of  Technology 
Sea  Grant  Program.^  ^ 
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Characteristics  of  Coastal  Stations  Located  1-3  Miles  (1.6-4.  8 kms) 
from  Shore  between  Cape  Ann,  Massachusetts  and  Rye,  New  Hampshire 

(Mulligan  and  deLara,  1974) 


Salinity  ( loo) 

Temperature  (°C) 

Depth  of  Euphotic  Zone  (m) 

3 

Chlorophyll  a (mg/m  ) 


9.  3-33.3 
1 - 18.  1 
5-20 
0.  3 - 15.  5 


Figure3-17.  Map  of  Northern  Massachusetts  and  New  Hampshire 

Coastal  Region 

A.  Rye,  New  Hampshire,  B.  Newburyport,  Massachusetts, 

C.  Annisquam,  Massachusetts,  D.  Salem,  Massachusetts, 

Stations  2-6  occupied  1973-1974. 


The  coastal  region  from  Rye,  New  Hampshire  to  Cape  Ann  is 

part  of  the  Gulf  of  Maine.  The  coastline  of  the  region  consists  mostly 

of  low,  sandy  barrier  beaches  which  protect  extensive  spartina  marshes 

and  estuarine  areas.  The  coastal  region  contains  the  Merrimack  River, 

the  major  river  between  Cape  Cod  and  Cape  Elizabeth  (Saco  River).  The 

Merrimack  River  drainage  basin  comprises  eight  percent  of  the  entire 

(23) 

drainage  basin  of  the  Gulf  of  Maine  and  dominates  the  physical,  chemical 
and  biological  characteristics  of  the  surrounding  coastal  region. 

(2.4) 

Average  annual  rainfall  in  the  region  is  43  inches  (1.09m)  which 
is  distributed  rather  uniformly  over  the  entire  year.  River  runoff  exhibits 
a definite  season  1 cycle  coinciding  with  the  period  of  major  snow  melt 
during  the  spring  (Table  3-2). 


TABLE  3-2 

Water  Discharge  Data  of  Merrimack  River  Below  Concord  River  at 
Lowell,  Mass,  in  CFS  (cubic  ft.  per  second) 


1. 


Date 

1945/46 

1 1970/7  i 2 

October 

4266 

3093 

November 

6047 

4257 

December 

10041 

4969 

January 

8856 

4169 

F ebruary 

7148 

5235 

March 

19200 

9843 

April 

9739 

21700 

May  ( 

11869 

16510 

June 

8536 

3511 

July 

2695 

1567 

August 

4181 

2109 

September 

2074 

erage  discharge  for 

48  year  period  ending  1971  = 

7055  CFS 

FORD,  W.  L.  1947. 

Hydrography  of  the  Western 

Atlantic.  No.  3 The 

distribution  of  the  Merrimack  River  effluent  into  Ipswich  Bay. 
Hole  Oceanographic  Institution  Technical  Report.  2 3p. 


W oods 


2. 


ANON,  1971.  Water  Resources  Data  for  Massachusetts,  New  Hampshire  , 
Rhode  Island  and  Vermont.  Part  I.  Surface  Water  Records.  Part  2. 
Water  Quality  Records.  U.  S.  Dept,  of  the  Interior.  Geol.  Survey.  401 


pp. 
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Winds  follow  the  same  general  patterns  as  elsewhere  in  the  Gulf  of 
Maine  (25)(26)  mostly  from  the  west,  with  considerable  seasonal  variation: 
in  Fall  and  Winter,  mainly  from  north  to  west;  in  Spring,  winds  are  variable 
and,  during  the  Summer,  winds  blow  from  the  sourthwest  and  southeast 
quadrants. 

Non-tidal  currents  tend  to  move  the  surface  coastal  waters  in  the 
region  counter-clockwise,  to  the  southwest  along  the  coast;  bottom  waters 
move  either  to  the  southwest  or  toward  the  shore  according  to  Bigelow 
Day  ^ ^ , and  Graham^  ^ . Graham  also  noted  that  bottom  water  which 

originated  in  the  northern  Gulf  of  Maine  moved  during  the  summer  through 
Jeffrey  and  Scantum  Basin  and  toward  the  shore  near  the  mouth  of  the 
Merrimack  River. 

Winds  can  generate  surface  currents  up  to  2%  of  wind  velocity,  and 
consequently,  non-tidal  currents  are  stronger  in  the  fall  and  winter  because  of 
higher  wind  velocities  from  the  northern  quadrant.  Currents  are  increased 
during  the  spring  due  to  increased  snow-melt  runoff  piling  up  along  the 
shoreline.  Winds  from  the  southern  quadrant,  common  in  the  late  spring  and 
summer,  reduce  and  even  reverse  the  slow-moving,  counter-clockwise,  surfac 
current  along  the  coast.  The  most  prominent  feature  of  the  water  circulation 
along  the  western  coast  of  the  Gulf  of  Maine  is  upwelling 

Mulligan  and  deLara^^^  have  recently  reviewed  the  Phytoplankton  and 
hydrographic  conditions  in  the  region  (Figure  3-17). 

Salinities  measured  from  August  1971  to  July  1972  ranged  from  9.  3 to 

33.3  0/00.  Highest  salinities  were  recorded  in  January  and  February. 

(?&) 

Manohar-Mararaj  and  Beardsley  estimated  that  the  Merrimack  River 
contributed  over  90%  of  the  fresh  water  in  Massachusetts  Bay.  Ice-flows 
and  low  salinity  surface  water  from  the  Merrimack  River  were  observed  six 
miles  offshore  and  could  be  seen  moving  southeastward  around  Cape  Ann  into 
Massachusetts  Bay  (Mulligan,  unpublished).  Consequently,  the  surface  waters 
off  and  to  the  south  of  the  Merrimack  River  displayed  the  lowest  salinities. 
During  the  springtime,  the  entire  water  column  off  Annisquam  had  reduced 
salinites  when  compared  with  conditions  present  at  coastal  stations  to  the  north. 

Water  temperatures  ranged  from  1-18.  1°C  over  this  period  and  were 
slightly  lower  in  the  northern  portion  of  the  region  throughout  the  year.  The 
coastal  waters  to  the  north  of  the  Merrimack  River  (to  a depth  of  20  meters) 
were  generally  well  mixed.  Thermal  stratification  developed  during  the 
warmer  months  and  became  well  established,  especially  at  the  mouth  of  the 
Merrimack  River  and  in  the  waters  to  the  south. 
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The  depth  of  the  euphotic  zone  varied  from  5-20  meters  with  the 
lowest  transparencies  occurring  during  the  spring  and  fall.  Silt  contributed 
by  the  Merrimack  River  was  largely  responsible  for  the  minimal  spring 
transparencies. 

Tables  3^3  , 3-4  , and  3-5,  are  reprinted  in  their  entirety  as 
taken  from  U.  S.  Coastal  Pilot  #1  NOAA  for  the  convenience  of  the  reader 
in  familarization  of  the  region,  climate,  mean  surface  water  temperature/ 
density  meterological  data. 

Figures  3-18  through  3-24  are  cuts  from  Nautical  Charts  13274,  13279 
and  13281  showing  current  R/V  berthing  location,  (13281)  cite  of  initial  apparatu 
shakedown  data  and  hydrographic  survey  information  of  the  region. 

Significance  of  the  bold  insets  are  as  below: 


Chart 


13281 

13279 

13279 

13279 

13274 

13274 

13281 


R/V  berth  near  Blynman  Canal 

R/V  berth  and  Annisquam  River 

Gloucester  Harbor,  Raw  sewage  outlet 

Deep  (56')  "clear"  water  hole  off  Coffins  Beach 

Deep  open  waters  off  Castle  Neck 

Merrimack  River  test  site. 

Sheltered  shakedown  area. 
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TABLE  3-3 


CLIMATOLOGICAL  TABLE 

DOSTON.  MASSACHUSETTS  <43*33'N. . Tl  OJ'W. > Elevation  IS  (L  (4.51m) 
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33.0 
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30 
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44.6 

56.3 
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39.3 

56.7 

30 
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31.3 

40.1 

50.  1 

39.  3 
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12 
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34 
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30 

IS 
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0.  63 

0.  35 

0.  96 

1.  72 

1.03 

0.  33 

23 

Maximum  In  24  hrs.  (Inches) 

2.07 

2.  ea 

4.  13 

2.31 

5.  74 

2.  46 

2.  42 

1.  40 
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Only  marine  radiobeacons  hove  been 
colibroled  for  surface  use.  Limitations  on 
the  use  of  certain  other  radio  signals  01 
aids  to  marine  navigation  can  be  found  in 
the  U.S.  Coast  Guard  Light  Lists  and  Oefenw 
Mopping  Agency  Hydrographic  Center 
Publication  117  (A&8). 

Rodio  direction-finder  bearings  to  com- 
mercial broadcasting  stotions  are  sobje 
error  and  should  be  used  with  caution. 

Stotion  positions  are  shown  thus:  J 
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3.  4 Q & A on  Fundamental  and  Planned  Experimental  Techniques 

During  the  course  of  the  Task  I performance  period  a number  of 
questions  were  directed  to  our  attention  regarding  both  fundamental  and 
planned  experimental  techniques  from  ONR  and  other  interested  parties. 

This  section  compiles  and  directly  addresses  these  questions  in 
dialog  format. 

Question  # 1 

What  effect  will  surface  waves  have  on  the  Raman  subsurface  temperature 
and  salinity  measurement  if  the  measurement  is  made  from  an  airborne  platform? 
Answer  #1 

The  major  surface  effects  relate  to  (a)  the  differential  transmission  of 

orthogonal  linear  polarizations,  (b)  the  loss  of  depth  resolution  due  to 

refraction  produced  by  transmission  through  the  air  -water  interface,  and 

(c)  the  lateral  errors  also  produced  by  refraction.  The  work  on  the  present 

contract  included  calculations  of  the  errors  produced  by  polarization  unbalance 

and  the  results  which  show  less  than  a 1.0%  effect,  which  is  equivalent  to  a 

1 . 0°C  temperature  error  if  uncorrected.  This  information  is  shown  in  Figure 

3-3  of  this  report.  The  depth  resolution  effect  and  the  lateral  error  were  also 

investigated  during  the  present  contract  and  the  results  are  shown  in  Figures 
3-5  and  3-6. 

Question  § 2 

Will  multiple  scattering  in  natural  waters  be  the  limiting  factor  in  the 
measurement? 

Answer  #2 

It  appears  that  multiple  scattering  will  not  be  the  limiting  factor  for  at 
least  the  first  15  to  20  attenuation  lengths.  This  is  because  (a)  the  very  sharply 
forward  peaked  volume  scattering  function  limits  the  beam  spread  to  the  order 
of  3 to  5°  and  (b)  small  angle  scattering,  i.  e.  5°  or  less,  exhibits  negligible 
depolarization.  For  this  assessment,  the  work  on  the  present  contract  drew 
heavily  upon  the  theoretical  work  of  A.  Gordon^  ^ and  the  experimental  results 
of  S.  Q.  Duntley^1  In  addition,  original  calculations  were  performed  by  CGC 
personnel  during  this  reporting  period  to  verify  that  Mie  scattering  at  small 
forward  angles  has  negligible  depolarization  and  also  to  calculate  the  scaling  laws 
that  determine  the  effect  of  the  fluctuation  in  depolarization  on  the  error  analysis 
of  the  measurement. 
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Question  #3 

To  what  extent  will  laser  power  fluctuations  limit  the  measurement? 

Answer  #3 

It  is  standard  practice  in  laser  Raman  measurements  to  monitor  the 

laser  power  and  normalize  the  Raman  signals  to  the  laser  power.  In  addition, 

real  time  computer  techniques  can  be  employed  to  repetitively  scan  spectral 

data  and  coherently  add  the  Raman  spectra,  with  noise  sources  such  as  laser 

power  fluctuations  being  effectively  subtracted  out.  Such  techniques  were 

(31) 

utilized  by  CGC  in  the  Raman  combustion  diagnostic  system  that  was 
installed  in  the  Air  Force  Aero  Propulsion  Laboratory  and  produced  Raman 
data  with  a precision  exceeding  1%.  Both  techniques  will  be  employed  in  the 
Raman  subsurface  temperature  and  salinity  experiments.  It  is  thus  not 
expected  that  laser  power  fluctuations  will  be  a limiting  factor  in  the  experiments. 

Question  #4 

In  the  orthogonal  polarization  ratio  technique  experiments  can  the  two 
photomultipliers  be  balanced  to  the  extent  required? 

Answer  #4 

Photon  counting  techniques  will  be  used,  and  in  such  a "0"  or  "1"  mode 
of  operation  the  constancy  of  the  photoefficiency  of  the  photocathode  surface  is  the 
important  factor,  not  the  gain  of  the  dynode  chain.  In  addition  the  polarization 
channels  can  be  optically  alternated  relative  to  the  two  photomultipliers  to  enable 
cross  calibrations  on  a regular  basis. 

Question  #5 

How  will  the  Raman  subsurface  temperature /salinity  experiment 
compensate  for  ambient  daylight  background  illumination,  especially  considering 
the  fact  that  ambient  upwelling  ocean  water  irradiance  is  usually  polarized? 
Answer  #5 

Operating  in  the  pulsed  laser  mode  the  experiment  will  provide  for 
sampling  of  the  ambient  upwelling  ocean  water  irradiance  background  during 
the  time  interval  between  laser  pulses.  With  two  polarization  channels  the 
ambient  background  in  the  orthogonal  polarizations  will  be  measured  and 
separately  applied  as  a background  correction  to  the  Raman  measurements. 
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Question  #6 

It  has  been  said  that  because  the  Raman  effect  is  a "weak"  interaction 
it  may  not  produce  sufficient  signal  for  the  remote  measurement  of  subsurface 
ocean  water  temperature  and  salinity.  Please  comment. 

Answer  #6 

Atmospheric  scientists  who  have  used  the  Raman  technique  to  measure 
atmospheric  water  vapor  profiles  are  familiar  with  the  fact  that  the  gas 
phase  vibrational  Raman  cross-section  is  usually  only  one  part  in  one  thousand 
or  less  of  the  on-frequency  Rayleigh  backscatter.  In  liquid  phase  Raman 
applications  not  only  is  the  molecular  density  several  factors  of  10  higher  than 
that  of  gas  phase  molecular  density,  but  also  the  per  molecule  Raman  cross- 
section  is  considerably  higher.  For  example,  the  Raman  vibrational  O-H  band 
backscatter  coefficient  from  distilled  water  is  essentially  equal  to  the  on- 
frequency  "Rayleigh"  backscatter  coefficient.  As  part  of  the  work  on  the  present 
contract  the  current  status  of  the  information  available  on  the  backscatte  ring 
coefficients  - both  Raman  and  on-frequency-was  researched.  The  results  are 
shown  in  Table  2-1  of  this  report. 

Question  if 7 

Will  salinity  gradients  interfere  with  the  measurement  of  temperature? 

Can  salinity  also  be  measured  by  using  the  Raman  spectrum?  How  accurate 
are  the  measurements? 

Answer  #7 

If  the  depolarization  ratio  is  spectrally  resolved  over  the  liquid  water 
O-H  band  vibrational  Raman  band,  it  is  observed  that  in  the  lower  frequency 
portion  of  the  band  corresponding  to  the  hydrogen  bonded  "ploymer-like" 
structures  the  depolarization  ratio  is  essentially  independent  of  salinity  and 
is  only  a function  of  temperature.  It  is  also  observed  that  in  the  higher 
frequency  portion  of  the  band  which  corresponds  to  non-hydrogen  bonded  or 
"monomer -like"  structures,  the  depolarization  ratio  is  a function  both  of 
salinity  and  temperature. 

Spectral  analysis  of  the  Raman  depolarization  ratio  thus  enables 
independent  determination  of  temperature  and  salinity.  The  sensitivity  of 
the  temperature  measurement  requires  a 1%  measurement  of  the  Raman 
depolarization  ratio  to  obtain  a 1°C  temperature  precision.  A three  part  per 
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thousand  (ppt)  accuracy  in  salinity  requires  0.8%  and  0.4%  Raman  measure- 
ment accuracies  in  the  salinity  ranges  of  0 to  20  ppt  and  20  to  40  ppt 
respectively.  It  should  be  noted  that  a temperature  change  of  1°C  and 
a salinity  change  of  3 ppt  each  has  the  same  relative  effect  on  sound  speed. 


Question  #8 

What  else  can  a high  repetition  rate  pulsed  laser  Raman  lidar  system 
measure,  if  it  is  fully  developed  to  measure  water  temperature  and  salinity? 
Answer  #8 

The  Raman  backscatter  from  the  water  can  also  be  analyzed  to  determine 
transmission  as  a function  of  depth.  In  addition  the  on-frequency  return  can 
be  utilized  as  a bathymeter  and  as  a wave  profiler.  Also,  if  employed  from 
an  airborne  platform  the  same  basic  system  can  measure  atmospheric 
temperature,  water  vapor,  transmission,  and  wind  profiles.  The  above 
described  multi-mode  utilization  involves  mainly  additional  are  for  data 

interpretation  with  relatively  minor  additional  hardware  costs 


Question  #9 

What  depth  resolution  is  possible? 

Answer  #9 

The  depth  resolution  depends  on  the  laser  pulse  duration,  the  receiver 
band  width  and  the  properties  of  the  water  transmission  medium,  especially 
the  refractive  properties  of  the  surface,  i.  e.  the  air-water  interface.  Off-the- 
shelf  laser  and  receiver  systems  have  demonstrated  30  cm  resolution  in 
bathymetric  lidar  systems,  which  are  generically  similar  to  a temperature / 
salinity  lidar  system.  The  effect  of  the  surface  was  discussed  in  the  Answer 
to  Question  #1,  and  the  results  obtained  from  the  investigation  during  this 
reporting  period  are  shown  in  Figure  3-5. 

Question  § 1 0 

Can  the  laser  Raman  method  detect  ice  and  measure  its  temperature? 
Answer  #10 

Yes.  The  laser  Raman  technique  has  been  shown^^  capable  of  detecting 
ice  and  measuring  ice  temperatures  in  the  range  from  0°C  to  - 50°C.  The 
physical  principle  is  the  measurement  of  Stokes  to  anti-Stokes  line  ratios  in 
ice,  which  are  functions  of  temperature. 
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Question  #11 

Can  a laser  Raman  system  be  built  to  produce  measurements  with  an 
accuracy  of  1%  or  better  in  a real  world  field  environment  as  contrasted  to 
what  might  be  achieved  in  a controlled  laboratory  environment? 

Answer  #11 

Laser  Raman  systems  have  been  built  for  remote  measurements  in  field 

environments  and  have  achieved  accuracies  under  such  conditions  of  1%  or 

(21 ) 

better.  Such  a system  is  the  laser  Raman  combustion  diagnostic  system 

that  was  installed  and  operated  by  CGC  in  the  Air  Force  Aero  Propulsion 

Laboratory.  The  combustion  diagnostic  system  has  a "measurement  index" 

- 1 4 

(i.  e.  ratio  of  laser  photons  transmitted  to  Raman  photon  collected)  of  ~10 
-12 

compared  to  ~10  for  a laser  Raman  water  temperature/salinity  measure- 
ment system.  On  the  basis  of  "measurement  index"  a water  temperature/ 
salinity  field  system  is  a more  conservative  undertaking  than  what  has  already 
been  demonstrated  and  operated  under  field  conditions. 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


This  section  summarizes  the  various  major  conclusions  and 
recommendations  that  were  developed  during  the  course  of  the  new  work 
performed  on  the  subject  contract.  The  following  items  are  arranged 
in  the  approximate  order  in  which  they  are  discussed  in  the  Technical 
Summary  Report. 

Liquid  Water  Raman  Cross-Section.  A new  value  for  the  Raman 
cross-section  for  liquid  water  appeared  in  the  recent  literature^ 
that  is  a factor  of  5 greater  than  that  previously  reported.  The 
implications  of  this  higher  value  should  be  factored  into  future  system 
design  considerations. 

Ice  Measurements.  It  has  also  been  reported  in  the  recent 
literature^  ^ that  the  temperature  of  ice  can  be  measured  with  Raman 
scattering.  A multi-purpose  Raman  instrument  is  conceivable  which 
remotely  measures  ice  temperature  in  addition  to  liquid  water  temperature 
and  salinity. 

Surface  Wave  Refraction  Effects.  It  was  concluded  that  since  surface 
wave  refraction  effects  are  of  a relatively  minor  importance  to  botli 
polarization  measurements  and  to  the  depth  and  lateral  resolution  that  is 
attainable,  no  attempt  should  be  made  in  the  initial  experiments  to  improve  the 
accuracy  of  the  measurements  by  correlation  on  a laser  pulse  to  pulse 
basis  of  the  corresponding  wave  angles.  Even  for  worst  case  Stokes 
limit  waves  of  H/L  = 0.  143  the  uncorrected  errors  are  0.6°C  temperature, 

2.8  parts  per  thousand  salinity  and  0.3%  in  depth  resolution. 

Depolarization  by  Transmission  in  Natural  Waters.  An  important 
conclusion  is  that  the  polarization  content  of  a beam  of  light  is  preserved 
essentially  until  the  last  photon  has  been  absorbed,  i.  e.  the  attenuation 
coefficient  is  very  much  larger,  the  order  of  100  times  larger,  than  the 
depolarization  coefficient  ^Experiments  should  be  carried  out  and  are  being 
planned  to  study  the  extent  of  the  natural  variability  of  the  proportionality  between 
P and  K.  This  information  will  define  the  residual  temperature  and  salinity 
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uncertainty  that  can  be  expected  from  a Raman  measurement.  Estimates 

are  that  with  a 10%  uncertainty  in  the  value  of  the  proportionality 

constant  between  P and  K an  uncertainty  of  0.  1 C in  temperature  and 

s 

0.7  parts  per  thousand  in  salinity  will  result.  Experiments  should  concentrate 
on  the  variability  and  the  uncertainty  in  the  ratio  of  P g to  K for  various 
types  of  natural  water. 

Depolarization  by  Transmission-Theory.  Basic  Mie  scatter  calculations 
confirmed  the  relatively  weak,  i.  e.  Pg  w 0.01  K depolarization  that  had 
been  experimentally  observed  in  the  forward  dir ectioJ.1 1 ^he  conclusion 
from  the  work  is  that  direct  experimental  measurements  of  depolarization 
in  various  types  of  sea  water  would  be  more  useful  than  further  complex 
and  idealized  calculations. 

Differential  Attenuation  and  the  Need  for  Multiwavelength  Polarization 
Detection.  A simultaneous  multiwavelength  detection  or  wavelength 
scanning  detection  capability  is  required  if  a reasonable  depth  capability 
is  to  be  obtained  even  using  a polarization  ratio  method.  An  experimental 
wavelength  scanning  capability  with  a wavelength  resolution  of  a few 
tenths  of  nanometers  is  necessary. 

Also,  there  is  a requirement  for  a tunable  laser  source  to  position 
the  Raman  water  band  at  or  near  the  minimum  of  the  diffuse  attenuation 
(K)  vs.  wavelength  curve.  Both  considerations  were  factored  into  the 
experimental  plan. 

Non-Raman  Photon  Noise  Sources.  The  ambient  background 
illumination  should  be  separately  sampled  between  laser  pulses  in  both 
polarizations.  Also,  experiments  will  be  designed  under  the  proposed  plan 
to  measure  the  laser  induced  fluorescence  to  Raman  ratio  as  a function  of 
of  laser  excitation  wavelength  for  a variety  of  coastal  waters. 

Experimental  Apparatus  Design.  The  conclusion  is  that  suitable 
laser  Raman  equipment  is  now  available  at  CCC  for  conducting  meaningful 
experimental  studies.  An  especially  important  aspect  of  the  experimental 
apparatus  design  is  the  availability  of  an  operational  research  vessel, 
the  Makai,  which  is  currently  outfitted  to  conduct  laser  Raman  experiments 
in  coastal  waters. 


Preliminary  Experimental  Raman  Data  Obtained  From  Natural 
Waters.  Temperature  measurements  from  Raman  spectra  were  obtained 
from  natural  waters  aboard  the  Makai  which  agreed  with  laboratory 
calibrations.  The  fluorescence  to  Raman  ratio  in  the  waters  selected 
for  experiments  was  very  low. 

Experimental  Field  Test  Plan.  The  coastal  areas  in  the  Gulf  of 
Maine  region  offer  a sufficient  variety  of  well  documented  coastal  waters 
so  that  meaningful  experiments  can  be  carried  out  in  this  region. 
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APPENDIX  A 


PHYSICAL  AND  MATHEMATICAL  JUSTIFICATION 


1.0  INTRODUCTION 

This  appendix  describes  the  physical  basis  and  the  mathematical  justification 
for  the  airborne  subsurface  temperature  area  profiles. 

2. 0 BACKGROUND 

Liquid  water  is  one  of  the  most  common  substances  on  the  earth,  it  is  necessary 
for  life  and  is  of  wide  practical  utility.  Because  of  its  basic  importance  to  man, 
water  has  been  the  subject  of  considerable  scientific  inquiry  and  its  spectroscopy 
is,  in  general,  well  understood. 

The  H^O  molecule  is  highly  polar  and  possesses  a large  dipole  moment.  This 
fact  together  with  its  atomic  composition  allows  the  H2O  molecule  in  the  liquid 
state  to  interact  with  other  nearby  H^O  species  and  form  hydrogen  bonds.  An 
individual  H2O  molecule  may  participate  in  a maximum  of  four  hydrogen  bonds  - 
one  for  each  H-atom  and  two  for  the  O-atom. 

The  fact  that  liquid  water  contains  hydrogen  bonding  is  apparent  from  its  infrared 
and  Raman  spectroscopy.  The  fundamental  vibration  frequencies  due  to  the  O-H 
bond  occur  at  3652  cm-*  (symmetric  stretch)  and  3750  cm"*  (asymmetric  stretch). 
The  O.  . .H  hydrogen  bond  that  occurs  in  the  liquid  state  perturbs  the  fundamental 
frequency,  so  that  the  infrared  and  Raman  spectra  due  to  O-H  in  liquid  water  show 
spectral  broadening  and  shifts  to  lower  frequencies  i>-  the  3400  to  3500  cm"  * range. 
Since  the  hydrogen  bond  is  relatively  weak,  having  an  energy  of  only  a few  kcal/ 
mole,  the  extent  of  hydrogen  bonding  in  liquid  water  is  very  sensitive  to  temper- 
ature. By  monitoring  the  Raman  spectrum  of  water  in  the  O-H  stretch  region, 
this  variation  in  structure  due  to  the  mobility  of  the  hydrogen  bonds  may  be 
observed  and  provides  a means  for  measuring  the  temperature. 

Walrafen  (References  Al,  A2,  A3)  was  one  of  the  first  to  observe  the  fact  that  the 
Raman  spectrum  of  liquid  water  in  the  O-H  stretch  region  shifted  toward  longer 
wavelengths  with  increasing  temperature.  It  was  also  observed  that  a much 
smaller  red  shift  occurs  as  ionic  species,  e.g.  Cl-ions,  are  added  to  the  water. 
Both  of  these  phenomena,  as  well  as  additional  spectroscopic  information  and  other 
data  have  been  used  to  support  the  hypothesis  that  liquid  water  contains  H2O 
clusters  held  together  by  substantial  hydrogen  bonding  in  rapid  equilibrium  with 
monomeric  H2O  molecules.  As  the  water  temperature  increases,  this  equilibrium 
shifts  to  a larger  fraction  of  unbonded  units.  Analogously,  as  the  ionic  concen- 
tration is  increased,  the  hydrogen  bonded  structures  would  be  decreased  since 
more  H2O  molecules  would  be  needed  to  cluster  around  the  ionic  species. 

3.0  TWO  COLOR  METHOD 

This  model  is  shown  in  a heuristic  manner  in  Fig.  A-l.  Liquid  water  is  considered 
as  being  composed  of  two  species,  polymers  which  are  hydrogen  bonded  and 
monomers  which  are  not,  both  in  rapid  equilibrium  with  each  other.  Each  species 
has  its  individual  Raman  spectrum  proportional  to  its  concentration  as  shown  in 
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A-l.  Raman  Spectrum  and  Structure  of  Liquid  Water. 


the  schematic;  the  total  observed  Raman  spectrum  being  composed  of  the  sum. 
Simultaneous  measurement  of  the  Raman  spectra  at  the  two  wavelengths  A ^ and 
7\,  ^ provides  a measure  of  the  concentration  ratio,  and  through  the  equilibrium 
constant,  the  temperature.  Such  a two-species  mixture  model  of  liquid  water 
has  previously  been  successfully  used  in  interpreting  sound  absorption  measure- 
ments (Reference  A4  ),  and  infrared  (References  A5  and  A6)  and  Raman  (Refer- 
ences A2  and  A7)  spectroscopic  studies. 

Chang  and  Young  (Reference  A8)  experimentally  demonstrated  the  feasibility 
of  the  two  color  technique  in  a series  of  laboratory  measurements.  They  used 
a N2  pumped  dye  laser  to  obtain  a laser  source  with  a wavelength  between  4200 
and  4600  A.  In  an  extensive  series  of  experiments  Chang  and  Young  studied 
the  sensitivity  of  the  temperature  measurement  to  the  location  of  the  two 
detected  colors  and  to  changes  in  water  salinity. 

Some  of  the  major  conclusions  of  this  feasibility  study  were:  (1)  over  a wide 
range  of  temperature  (0-50°C)  the  technique  was  linear  with  temperature,  (2) 
two  colors  could  be  found  which  maximized  the  temperature  sensitivity  and 
minimized  the  salinity  effect,  (3)  a measurement  sensitivity  of  1%  per  C°  was 
demonstrated  and  (4)  a two  color  field  system  for  water  temperature  measure- 
ment was  found  to  be  feasible.  The  differential  attenuation  of  the  Raman  band 
at  depths  corresponding  to  several  attenuation  lengths  was  identified  as 
requiring  an  independent  measurement  of  transmission  as  a function  of  wavelength 
as  a calibration  requirement  in  the  two  color  method. 

4.  0 DEPOLARIZATION  METHOD 

As  mentioned  previously,  liquid  water  may  be  considered  as  a mixture  of 
species  which  are  completely  hydrogen  bonded  and  those  which  are  not.  The 
former  possess  symmetric  structure  and  vibrations  and  would  be  expected  to 
preferentially  scatter  polarized  radiation.  Hence,  in  the  same  way  that  liquid 
water  structure  is  reflected  in  the  shape  of  the  Raman  spectrum,  similar 
information  is  contained  in  the  polarization  of  the  scattered  radiation. 

Porto  (Reference  A9)  has  shown  that  the  total  Raman  intensity,  I,  as  a function 
of  scattering  angle,  6 , is  given  by  the  expression 

1 ( 0)  = — 4T"  C45  a 2 + 7 ^ 2] 


where  I is  the  intensity  of  the  incident  radiation  and  a and  0 are  the  isotropic 
and  unisotropic  part  of  the  differential  polarizability  tensor.  a and  0 are 
related  to  the  molecular  polarizability  and  are  temperature  dependent.  For 
■bservation  of  back  scattered  radiation,  0=  0.  For  linear  polarized  incident 


« jation,  a depolarization  ratio 


can  be  defined 


3 0 


45  a 


4 0 
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wnere  I,  and  Itl  are  the  parallel  and  perpendicular  polarized  polarization 
intensities.  As  has  been  shown  by  Murphy  and  Bernstein  (Reference  A7) 
the  bulk  of  the  Raman  return  for  water  is  composed  of  IM.  Hence />  is 
a small  number  and  the  accuracy  of  such  a measurement  is  inherently  low. 

However,  circularly  polarized  light  can  be  used  for  Raman  scattering,  and  in 
an  analogous  manner  to  the  plane  polarized  radiation,  the  scattered  return 
will  contain  components  of  both  circular  polarizations.  A reversal  coefficient 
or  circular  depolarization  ratio,  />  , can  be  defined  and  related  to  the  linear 
depolarization  ratio  (Reference  A13j: 


P c 


when  1 x and  are  the  left  and  right  circularly  polarized  return  respectively. 
Since^  <3/4,  /®c  <6  and  is  a more  accurate  measurement  of  water 


temperature  than  p 


Chang  and  Young  (Reference  A8)  were  able  to  determine  the  circular  depolari- 
zation ratio  of  liquid  water  as  a function  of  temperature.  They  found  that  the 
ratio  was  approximately  linear  over  a temperature  range  of  about  10-45°C  and 
determined  an  average  sensitivity  of  about  1 %/C°  which  is  comparable  to  the 
maximum  sensitivity  observed  for  the  two  color  method.  In  addition,  they 
found  that  the  circular  depolarization  ratio  varied  with  salinity,  but  were  able 
to  find  wavelength  regions  where  the  salinity  effect  was  minimal  and  could  be 
neglected  for  all  but  the  most  accurate  measurements. 


The  performance  of  the  temperature  measurement  system  utilizing  the  depolari- 
zation ratio  will  be  degraded  by  any  xncalib rated  effect  which  depolarizes  either 
the  laser  beam  or  the  Raman  scattered  beams.  Particulates  in  the  ocean  will 
cause  some  depolarization  by  single  and  multiple  off  axis  scattering.  Single 
scattered  light  scattered  at  very  small  angles,  e.  g.  10"^  radians,  will  be 
collected  by  the  receiver.  Scattered  light  at  larger  angles  can  be  collected 
when  multiple  scattering  occurs.  The  off  axis  scattered  light  from  a particle 
in  general  will  be  elliptically  polarized.  Because  the  nature  of  the  scattering 
in  the  ocean  cannot  be  universally  characterized,  the  temperature  sensitivity 
of  the  system  will  be  degraded.  Chang  and  Young  performed  experiments  to 
determine  the  degradation  of  system  performance.  Their  estimate  was  that 
small ^ngle .scattering  would  introduce  an  absolute  temperature  uncertainty 
of  0.  3 C per  attenuation  length.  This  uncertainty  is  independent  of  laser 
power  or  integration  time. 

5.  0 NUMERICAL  CALCULATION 

The  number  of  Raman  photoelectrons  collected  per  second  by  an  optical  detector 
in  a pulsed  laser  back  scattering  system  can  be  expressed  by  the  following 
equation: 
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NPE  “ NLASER  NSCAT  ffRAMAN  AR^  * pe  *op  TXl  T^2 


where: 


N 

N 

N 


PE 

LASER 

SCAT 

RAMAN 

A R 
Q 


t 

pe 


c 

op 


T>  1 T ^2 


number  of  Raman  photoelectrons  detected 
- number  of  laser  photons  per  second 


density  of  molecular  scatters  of  a given  species 


Raman  scattering  cross-section  per  particle 
per  steradian 

range  resolution 

detection  solid  angle 

photocathode  photoelectric  efficiency 


optical  system  efficiency 

two-way  transmission,  T ^ j at  laser  wavelength 
T ^ ^ wavelength 


It  can  be  seen  from  the  above  that  if  the  transmitter  and  receiver  system 
parameters  are  known  and  if  the  Raman  scatterer  has  a known  cross-section 
and  density,  such  as  liquid  water,  then  a measure  of  the  Raman  scattering 
from  a single  range  cell  can  be  used  to  determine  the  two-way  transmission 
to  that  range  cell.  The  system  thus  acts  as  a single-ended  transmissometer . 

The  Raman  technique  for  transmission  measurements  in  the  atmosphere  has 
already  been  demonstrated  by  Leonard  and  Caputo  (Reference  A10). 

A calculation  was  carried  out  for  the  baseline  Raman  system  using  the  above 
range  equation  with  system  parameters  that  assume  state-of-the-art  well 
engineered  components.  The  parameters  assumed  are  listed  in  Table  A- 1 
with  a brief  rationale  for  the  particular  choice. 

The  product  of  the  system  parameters  shown  in  Table  A-l  yields  a total  Raman 
signal  return  of  3.5  x 10^  photoelectrons  per  second  with  the  assumption  of 
unity  two-way  transmission  loss.  As  has  been  previously  shown,  (Reference  All) 
in  a well  engineered  Raman  system  using  a pulsed  laser,  the  limiting  noise  can 
be  reduced  to  the  "shot  noise"  in  the  laser  induced  signal  itself.  Therefore,  the 
S/N  ratio  can  be  given  by  the  square  root  of  the  total  number  of  photoelectrons 
collected  in  a given  measurement  interval.  If  the  signal  in  the  above  calculation 
is  assumed  to  be  equally  divided  between  the  two  polarizations  a basic  S/N  of 
(1.75  x 10^)  1/2  or  418  would  be  available.  Based  on  a measurement 
sensitivity  of  1%  per  C°  (Reference  A8).  This  would  mean  that  a temperature 
precision  of  (1  /418)  x 100  « 0.3  C°  would  have  been  achieved. 

The  ability  to  achieve  depth  penetration  depends  critically  on  the  attenuation 
coefficient  of  the  particular  water  and  is  highly  variable  among  open  ocean. 


A5 


continental  shelf  and  coastal  waters.  Table  A-2  lists  the  various  types  of 
ocean  waters  with  values  of  attenuation  coefficients,  optimum  transmission 
window  ranges  and  the  rate  of  laser  beam  power  loss  in  dB  per  meter. 


TABLE  A- 1 

BASELINE  SYSTEM  PARAMETERS 


System 

Parameter 


^Laser 


N 


Scat 


a 

Raman 


Value 

Assumed 


Rationale 


2.3  x 10 


17 


photons 

sec 


-5  -1 

10  s te  r 


m 


-1 


This  value  corresponds  to 
0.  1 watt  of  average  power 
at  4600  K.  C ommercially 
available  N2 -pumped  dye 
lasers  have  this  capability. 


This  is  the  latest  published 
value  for  the  experimentally 
measured  total  Raman 
liquid  water  backscatter 
coefficient . (Refe rence  A8) 


A* 


£1 


1 meter 


ster 


Corresponds  to  minimum 
depth  resolution  set  by 
pulse  duration  of  laser. 


This  collection  solid  angle 
corresponds  to  a 30  cm 
effective  diameter 
collector  on  a platform 
300  meters  above  the 
scattering  volume. 


op 


« 

pe 


5 x 10"2 


Typical  overall  optical 
and  photoelectric 
quantum  efficiency. 


UNDERWATER  WINDOWS  FOR  VARIOUS  OCEAN  WATERS 
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SWDCO  (SURFACE  WAVF  DEPTH  CORRECTION) 


i 


APPENDIX  C, 

I 

Computer  program  to  determine  the  effect  of!  ocean  surface 
waveshape  on  Refractive  depth  penetration  correction 
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OEPTH  CORRECTION  SUMMARY  SHEET  FOP  A SINOSQIOAL  OCEAN  WAVE 
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DEPTH  CORRECTION  SUMMARY  SHEET  FOR  A TROCHO  T DAL  OCFAN  WAVE 


r 


r\ 

Vi' 


TH 
FBO» 


rH 


Z 

o 

i* 

E 

K 

w 


L 

z 

ft 

I 

t 

© 


l(u 

iO 

k 

a 

ft 

< 


It- 

z 

IL 

u 

tt 

in 

It 

E 


I 

z 

IM 

Z 


\ 

X 


m ©min  in  min  in  <t  <e  <j  r ri  mr)  mnrir!  m m w;w  w w w w w w njm 

I I ; I II  I'll  I I II  III  I I III  I I II  I I till  ij  I 

w — ACirww©«f>c  — crc  — wmoerxmmoee  w<ox  « w 
m s w o © n © m © c <t  n r, s ® w in  o <f  <r  w,«c  w w m ec  w m © © 
NNffxCoetivNc-orffcrwM'iM'iNSM'i-iii-xa) 

> e ® mjw  n|< 


**  w a -•  «*  m -a  h-  — w -a  — w i 

! I 


*o  <cjm ir  in  tr if.  if  <e  a c r>  m m m m r, r r.  m m w w w w w w cy  m m m 

oo  ©CO©  00000  O©  ©OOO  O©  O OO  O©  © ©©  © O © © 

I I jl  I I i l I i i I I l i l l I l I l i I |l  l i i l l l I l 

« <0N  N * - W < IT  OV  (O  >■  <i  O iTNriMO  (?  r;r,  IT,  « « cc  (l(y 
<feac-eooror,co.^r:c-OKinr  n © cc  w -c  r ffC'j 

— 'CNoeemm— ^cr-*ewp.r-e-.mcr-»wjinmw-«m~trin® 

• •I  ••••••••••••[  ••••  «•••!•••  • ••••• 

— — cvj-c-«rjr’<^NNNNCk«^  — rrr,'  — — co  in  r: 

1 ! ' L ' [ ; I 

hi£i4'fycincirir<»<}<>«rir;nrrDrrirr(yM,'r)rnr,n 
i ii  i i i ! t i i ii  i i i j i i i i i i t ii  t i i i i i n 


©©©mwor-m— r^«fOwwowwm<fcm.~eroNio  — mi"^* 
in  ©!©  w tr  © **  © >c  ® <c  m o in  © <c  m o>  c — — © -« w oc  ©me  ©a  m 
-•  mi©  in  a in  w © -*m  © in  © if  w e ~ — w m © -■  m m — o;  m mr,  — in 

>0  ~ win  O'  — w n >C  -*'win  t — w m mr  m r © trim  — K in  © m w ••< 


cc  © cc.  n r-  m m m c m m in  in  m m if  m in  m m © © © © © © © © m m m 
o ov © e c-  o © c ©jo  © © « oeccooccoo©©©©©©© 

i l |i  i i i !i  i i i i i i i ,i  i i i i i i i i i i t i t i i l 

e ce'e  m a c •©  <r  c-  © c: tr  r —in  © n n n m © -c/h-  in  c w!<r  mi  tr . © 
in  ® © <;  • er  n ® in  v.  — i u — © c n m © o c ® ,©  ® © m,m 

~ or  jN  C u;  m w m m >c cj  cl  r n o — w © cy  © c U O'  w o w o a <rj— 

w © ® ~ m m«c  -« n o c - — — — — — — © ini© 

[ j i [ 

r-  r-  n n >c  m m 'll  <c  w in  in  in  in  m in  in  in  in  in  m in  in  m m 
! i } i I j i i I 1 ! ] i j 1 i l I 1 1 1 1 1 • 1 1 i i i i 


<r\Cm^mmNcc<r<r©nc«'C®sm>cr-©©©ir 
® r--  m w <v  m in  m in  m <o  c © © r-  eo  o © a,  m m © a 

co,tooonN-  m m © ® © -•  — — — — m^ii'bjjj  mmcciIh 


oorcco^N © ee-w©  me  — — — — — — 

l! 


U wrinmr-rmp 


oooooo.oecocc'<  rcoooyoocoocotriinnin 
w n © <c  or  © w m o o © c © o ’©  o w m m © © © !©  o c o n w w m m 
o c o c © — — — cv  ir  © «c  ® o w a <j  * © netpeo  © c m © — !© 
© ©pOOCO  t'COl>OC«r<-'-----'Nn<l  m - ••  — — o 
*1  • * 1 ■ 


© o©o  © ©;©  o©  ooC'  oc  oc  © o©  o © o 

I 


© o © o o © o o ;o  o o in  in  o o © 
|iceoinoooooo-MtiMfioo 
io  o — — w m * ink  n n n f-  n ® © 
poooproopooooooo1 


K>©  © © O©  O O;© 


o mo  o o o 

j -- 

mm  w in 

V oo 

oooo 

< oo 

oooo 

• • 

• • • • 

oo 

oooo 

pooo 


oooo 


paooooootoooo 


o o o op  o o o 
oco oo  o o o 
o m o o.o  o o o 
w — cy  pi  U m >0  ® 


oooo 
C8 


o c 


I 

o 


-•  I . nc  . - - 


APPENDIX  D 


DEPOLARIZATION  CALCULATIONS 


This  appendix  summarizes  calculations  which  show  the  scaling 
of  the  depolarization  produced  by  small  angle  scattering  during 
transmission  of  light  through  natural  ocean  waters  as  a noise 
source  in  addition  to  the  depolarization  produced  by  the  water 
Raman  scattering. 

Consider  a beam  of  light  with  polarized  component  Ip^  and 
unpolarized  component  IQ  incident  upon  a transmission  medium 
characterized  by  a differential  depolarization  coefficient  Pf  defined 
by  the  following  expression: 


dl 


At  some  transmission  distance,  Xj,  the  polarized  component 
will  become 


I e-p*Xl 
Po 


and  the  unpolarized  component  will  increase  to  become 


I 

U1 


I + I (1  - e”  pi X 1 ) 
uo  P«> 


(The  above  argument  neglects  attenuation  and  assumes  a homogeneous 
medium. ) 

If  at  the  transmission  distance,  Xj,  the  Raman  backscattering  process 
is  now  allowed  to  occur  the  polarized  component  Ip  will  be  converted 
via  the  Raman  process  into  two  orthogonal  polarization  components,  1^ 
and  IB 

where  .PR  = IA  / IB  and  IA  + IB  = I 


The  resultant  polarized  and  unpolarized  components  present 
following  the  Raman  process  are  therefore 


1 -P 


R 


i+pR  pi 


= i.  + 


2 PR 


U1  1 + 'r  pl 


If  the  receiver  is  co-located  with  the  transmitter,  the  Raman 
scattered  light  undergoes  further  depolarization  during  transit 
through  the  water  from  the  transmission  distance,  X.,  back  to  the 
transmitter /receiver  location.  The  return  signal  is  thus 


I = I e-'*Xl 

Ps  P2 


I = I + I (1  - e"P*Xl) 

Uj  u2  p2 


Or  expressing  the  return  signal  in  terms  of  the  initial  conditions 


I = — I 
M2  2 


where  the  overall  depolarization  ratio  for  the  two-way  water  transmission 
and  the  Raman  process  may  be  states  as 


L 


appendix  e 

{This  research  note  will  appear  in  Geophysical  Research  Letters,  July  1977.) 


EXPERIMENTAL  REMOTE  SENSING  OF  SUBSURFACE  TEMPERATURE 
IN  NATURAL  OCEAN  WATER 


D.  A.  Leonard,  B.  Caput o and  R.  L.  Johnson 
Computer  Genetics  Corporation,  LaVeside  Office  Park,  Wakefield,  Massachusetts  01880 

F.  E.  Hoge 

NASA  Wallops  Flight  Center,  Wallops  Island,  Virginia  13337 


Abstract.  A remote  sensing  laser  Raman 
backscattering  technique  has  made  the  first 
Held  measurements,  from  a boat,  of  subsurface 
temperature  in  natural  ocean  waters.  These 
results  and  the  measurement  systems  rapidly 
developing  from  them  open  wide  areas  of 
exploitation.  The  teclyiique  utilizes  the  obser- 


vation of  the  3400  cm"  O-H  stretching  Raman 


band  of  liquid  water.  The  spectral  data  were 
analyzed  to  obtain  the  subsurface  water  temper- 
ature in  accordance  with  previous  laboratory 
studies.  The  importance  of  such  a measure- 
ment by  the  Raman  scattering  technique  is  that 
it  now  enables  depth  resolved  profiles  of  ocean 
water  temperature  to  be  obtained  by  remote 
sensing  from  airborne  platforms  that  can  be 
adapted  for  application  to  many  problem  areas 
of  physical  oceanography,  meteorology,  clima- 
tology, ocean  technology,  and  thermal  dis- 
charge management.  An  accuracy  of  j^Z°C  has 
been  obtained  in  the  initial  field  experiments. 

A more  refined  system  has  been  designed  which 
will  permit  accuracies  of  better  than  + l°C  in 
the  first  one  or  two  optical  attenuation  lengths 
of  depth. 


The  purpose  of  this  communication  is  to  re- 
port the  first  successful  depth  resolved  remote 
sensing  measurements  of  subsurface  ocean 
water  temperature  in  a boatbome  field  experi- 
ment. Such  measurements  have  strong  impli- 
cations for  research  and  technology  in  the  gen- 
eral areas  of  ocean  and  coastal  physics,  ocean- 
ography and  weather  and  climate  prediction  and 
in  thermal  pollution  management.  The  temper- 
ature measurement  was  obtained  by  spectral 
analysis  of  the  3400  cm*1  O-H  stretching  Raman 
band  of  liquid  water  in  accordance  with  previ- 
ously reported  laboratory  experiments 
[Walrafen,  1961:  Schwiesow,  1971;  Chang,  et  al, 

ittci: — 


The  physical  basis  of  the  measurement  is  the 
fact  that  (a)  liquid  water  exists  in  at  least  two 
forms,  monomers  and  dimers;  (b)  the  two 
forms  arc  in  chemical  equilibrium  as  a function 
of  temperature;  (c)  the  O-H  Raman  stretching 
frequency  is  significantly  different  for  the  mono- 
mer and  dimer  forms;  (d)  the  relative  concen- 
tration of  monomer  vs.  dimer  can  be  determined 
from  the  Raman  spectrum  and,  thus,  the  temp- 
erature can  be  inferred. 

The  major  advantage  of  the  Raman  scattering 
remote  sensing  technique,  especially  when 
adapted  to  airborne  versions,  1*  that  it  can 
rapidly  and  economically  measure  subsurface 
water  temperature  profiles  over  large  areas. 

As  such  it  has  direct  application  to  the  elucida- 
tion of  such  specific  temperature  related  prob- 
lems (among  others)  as  input  to  sea-fir  inter- 
action modeling,  the  physics  of  oceanic  meio- 
scale  systems,  behavior  of  ocean  fronts  and 
other  mixed  layer  dynamics,  synoptic  behavior 
of  sound  propagation,  monitoring  of  thermal 
discharges  and  validadon  of  thermal  discharge 
models,  etc. 


El 


A laser  Raman  lidar  system  was  assembled 
and  operated  to  obtain  Raman  spectra  of  ocean 
water  aboard  the  Computer  Genetics  Corpora- 
tion research  vessel  Makai.  The  purpose  of 
the  effort  was  to  demonstrate  remote  subsurface 
temperature  measurements  in  natural  waters. 
Although  laboratory  experiments  had  been  per- 
formed previously,  Raman  temperature  mea- 
surements of  natural  ocean  water  had  never 
before  been  obtained  under  in  situ  conditions 
and  there  were  unanswered  questions  as  to  the 
extent  to  which  naturally  occurring  inter- 
ferences, such  as  fluorescence,  would  place  a 
limit  on  the  accuracy  and  the  resolution  of  the 
laser  Raman  technique  especially  in  coastal 
zones.  An  experimental  program  has  been 
designed  to  answer  these  questions.  This  paper 
reports  initial  results. 

The  block  diagram  of  the  laser  Raman  system 
that  was  operated  aboard  the  Makai  is  shown  in 
Figure  1.  The  laser  source  was  a pulsed  nitro- 
gen laser  operating  at  a wavelength  of  337.1 
nanometers  and  producing  pulses  of  100  kilo- 
watts peak  power  with  an  effective  pulse  dura- 
tion of  10  nanoseconds  at  a pulse  repetition  rate 
of  500  Hertz. 


The  output  from  the  laser  is  passed  through  an 
interference  filter  which  pastes  the  337.  1 am 
liter  line  with  high  efficiency  but  blocks  spon- 
taneous emission  occurring  in  the  laser  gas 
discharge  m the  O-H  Raman  spectral  region 
near  3800  A.  An  attenuator  is  also  placed  in  the 
laser  beam  to  provide  a beam  power  attenuation 
of  the  order  of  X100  so  that  a photon  counting 
detection  system  previously  employed  in  atmos- 
pheric Raman  measurement  applications  system 
could  be  used  in  the  liquid  water  studies. 

As  shown  in  the  block  diagram  a bistatic 
arrangement  of  the  transmitter  and  receiver 
optics  functionr  to  define  a volume  of  subsurface 
water  from  which  the  laser  Raman  scattering  is 
detected.  Operation  in  such  a bistatic  mode  is 
adequate  for  experimental  studies  in  the  first 
few  tens  of  meters  of  depth.  This  configuration 
is  being  modified  for  deeper  penetration  applica- 
tion by  utilizing  electronic  range  gating  tech- 
niques and  a more  suitable  blue-green  laser 
wavelength. 

The  photons  collected  by  the  receiver  optics 
are  passed  through  a double  1/4  meter  focal 
length  scanning  spectrometer  having  a 0.  S nano- 
meter spectral  resolution.  The  wavelength 
scanning  of  the  spectrometer  is  controlled  by  the 
computer.  The  function  of  the  spectrometer  is 
to  produce  a spectral  scan  of  the  liquid  water 
Raman  O-H  stretching  band  with  sufficient 
resolution  so  that  the  monomer  and  dimer  com- 
ponents can  be  resolved  and  a temperature 
dete  r mined. 

The  Raman  photon  output  from  the  spectrom- 
■eter  is  detected  by  a RCA  1P128  photomulti- 
plier and  amplified  by  an  Ortec  Model  454 
timing  filter  amplifier.  The  conditioned  Raman 
photon  signal  and  a gate  pulse  suitably  synchron- 
ised and  delayed  with  respect  to  the  laser 
firing  are  combined  in  a coincidence  circuit  the 
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output  from  which  is  recorded  using  direct 
memory  sects*  (DMA)  techniques  with  a Data 
• General  NOVA  1200  computer.  The  number  of 
pulses  occurring  during  the  same  data  taking 
interval  is  also  recorded.  The  ratio  of  the 
Ratnan  photons  to  laser  pulses  recorded  over  a 
given  time  interval  is  the  basic  data  obtained. 

The  data  taking  is  fully  automated.  The 
computer  can  be  directed  through  the  teletype 
to  specify  the  number  of  laser  pulses 
over  which  to  record  data  at  each  wavelength 
and  the  maximum  and  minimum  wavelengths  and 
the  wavelength  increments  of  each  spectral  scan. 
The  data  obtained  are  printed  on  the  TTY  in 
tabular  form  listing  "wavelength",  "data", 
"cumulative  data"  and  "drift".  Successive 
rapid  spectral  scans  to  construct  a "cumulative 
data"  record  allows  the  spectral  Information  to 
be  added  coherently,  thus,  Increasing  the  signal 
to  noise  ratio.  The  "drift"  data  record  is  a real 
time  diagnostic  that  allows  the  operator  to 
easily  recognize  during  an  experiment  the  occur- 
rence of  large  error  signals  as  soon  as  they 
may  occur. 

Initial  Raman  spectral  data  have  been  obtain- 
ed from  the  research  vessel  Makai  on  Novem- 
ber 3,  5 and  6,  1976  and  March  26,  27  and  30, 
1977  at  various  depths  from  the  surface  to  10 
meters  below  the  surface  in  the  Annisquam 
River,  a tidal  estuary  in  Massachusetts.  A 
typical  sample  of  such  data  obtained  at  a depth 
of  1 meter  below  the  surface  is  shown  in  Figure 
2.  The  measurements  were  made  under  a high 
tidal  flow  condition  in  waters  polluted  by 
industrial  and  ocher  urban  runoff.  It  was  esti- 
mated that  the  diffuse  absorption  coefficient 
was  of  the  order  K=0.3  m-1".  The  range  limi- 
tation in  this  prototype  bistatic  experiment  was 
a geometrical  one  due  to  the  inability  to  spa- 
tially resolve  the  Intersection  of  the  focussed 
spectrometer  slit  on  the  laser  beam  at  depths 
greater  than  10  meters. 

Figure  2 is  a plot  of  the  number  of  Raman 
photoelectrons  collected  per  0.5  nanometer 
spectral  interval  per  10,  000  laser  pulses  as  a 
function  of  wavelength.  Shown  are  two  sets  of 
data,  the  dashed  curve  obtained  as  an  instru- 
ment calibration  in  the  laboratory  with  a water 
sample  at  22°C  and  the  solid  curve  obtained  in 
the  Annisquam  River  water  at  6°C,  according 
to  "ground  truth"  temperature  readings  obtain- 
ed with  a pair  of  calibrated  mercury  thermom- 
eters. 

A temperature  can  be  obtained  from  the 
Raman  spectra  by  comparing  the  height  of  the 
band  contour  at  the  peak,  which  corresponds  to 
mostly  monomer  contributions,  to  the  height  of 
the  contour  at  the  inflection  point  on  the  short 
wavelength  side  which  corresponds  to  mostly 
dimer  contributions.  These  spectral  features 
are  labeled  monomer  and  dimer  respectively  on 
Figure  2.  - A scale  showing  the  height  of  the 
dimer  inflection  point  as  a function  of  temper- 
ature has  been  calculated  from  previously 
published  laboratory  data  [Walrafen.  1961]  and 
is  shown  on  Figure  2 as  a reference.  Bated  on 
this  scale  an  Annisquam  River  temperature  of 
4°C  is  inferred  from  the  spectra  and  is  con- 
sistent with  the  ground  truth  temperature  of 
6°C  to  within  the  shot  noise  limited  statistical 
precision  of  the  Raman  measurement  which 
was  + 2°C. 


* ol  initial  experimen- 

£nfl!.ld  ements,  the  performance  of  a 

fully  developed  airborne  User  Raman  water 
temperature  measurement  system  can  now  be 
estimated  with  a high  degree  of  confidence 
using  standard  lidar  engineering  calculations. 

f *“8  mi;P  °£  FigUre  3 b‘*°  calcu- 

lated  for  a typical  open  ocean  measurement 
situation  with  a mixed  layer  depth  of  100  meters. 
The  basic  parameters  on  the  map  are  the  num- 
ber  of  Raman  photoelectrons  collected  as  the 
ordinate  with  the  depth  as  measured  in  attenua- 
tion length*  displayed  on  the  abscissa.  The 

■ °£  precision  in 

c and  the  depth  in  meter*  a*  a function  of 
water  quality  are  also  shown  as  alternate  ordin- 
ate and  abscissa  quantities  respectively.  An 
fSef101*011117  mea3ured  fSlushe r and  Derr. 

1975J  cross-section  of  4.  5 ~x  10-29  cm*  7^51.- 
cule-sr  was  used  in  the  calculations. 

The  operating  lines  on  the  map  are  drawn  for 
various  values  of  the  parameter  M where 


M 


■(*)  (■*) 


(1) 


and  PL  1,  the  averaged  laser  power  in  watts, 

* b . depth  resolution  in  meters,  R is  the 
altitude  of  the  aircraft  in  kilometers  and  T is 
the  integration  time  of  the  measurement  in 
1 Jbus'  example,  for  a system 
with  Pl  = 0.2  watts,  Ah  * 10  meters,  R ■ 1 
Kilometer  and  T * 1 second,  the  map  parameter 
would  be  M = 0,01.  This  means  that  with  such 
a system  the  allowed  region  of  the  operating 
map  of  Figure  3 is  to  the  1 eft  of  the  line  labeled 
M = 0.  01.  A typical  operating  point  la  shown 
corresponding  to  a temperature  accuracy  of 
1 C and  a depth  penetration  of  60  meters. 

Also  shown  on  the  operating  map  Is  an  upper 
boundary  labeled  "volume  depolarization  inter- 
f*Jeuce  limit".  This  represents  the  most  sig- 
nificant interference  to  the  measurement  is 
caused  by  depolarization  due  to  transmission  of 
light  through  natural  water*.  Analysis  of  pub- 
lished data  [Duntley.  1971]  has  shown  that  the 
depolarization  effect  will  produce  uncertainty 
in  the  measurement  of  0. 1°C  per  attanustion 
length. 
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Fig.  2.  Raman  spectral  data  for  liquid  water 


Tig.  3.  Typical  system  operating  map  (100  meter  mixed  layer). 
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